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Abstract 

These notes are concerned with harmonic and holomorphic functions 
on Euchdean spaces, where "holomorphic" refers to ordinary complex 
analysis in dimension 2 and generalizations using quaternions and Clifford 
algebras in higher dimensions. Among the principal themes are weak 
solutions, the mean-value property, and subharmonicity. 
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1 Euclidean spaces 

1.1 Real and complex numbers 

Let R denote the real numbers, equipped with the usual arithmetic operations 
and ordering. We write Z and Z+ for the integers and the positive integers, 
respectively. The complex numbers are denoted C, and of course a complex 
number c can be written as a + ib, where a, b are real numbers and i"^ = —1. 

Recall that a real number b is said to be an upper bound of a subset A of R 
if a < 6 for all a € A. If A is a set of real numbers and c is a real number, then 
c is said to be the least upper bound or supremum of ^ if c is an upper bound 
for A and if c < 6 for every real number b which is an upper bound for A. It is 
clear from the definition that the supremum of A is unique if it exists. 

Similarly, a real number b' is said to be a lower bound for a subset A of R 
if 6' < a for every element a of A. If A is a subset of R and c' is a real number, 
then c' is said to be the greatest lower bound or infimum if b' < c' for every 
real number b' which is a lower bound for A. Again, it is easy to see from the 
definition that the infimum is unique if it exists. 

The completeness property of the real numbers, with respect to the ordering 
on the real numbers, states that every nonempty subset of R which has an upper 
bound also has a least upper bound. This is equivalent to saying that every 
nonempty subset of R which has a lower bound also has a greatest lower bound. 
The equivalence between the two statements can be shown using multiplication 
by —1, or by taking the supremum of the set of lower bounds of a subset of R 
to get a greatest lower bound and the infimum of the set of upper bounds of a 
subset of R to get a least upper bound, as appriopriate. 

The supremum and infimum of a subset A of R are denoted 



respectively, assuming that they exist. If f{x) is a real- valued function on some 
nonempty set E, then we may write 



for the supremum and infimum of f{x) on E, assuming that they exist. More 
precisely, these are the same as the supremum and infimum of the set 



(1.1) 



sup A, inf A, 



(1.2) 



sup /(a;), inf /(x) 



(1.3) 



{f{x) -.x&E) 



of values of /. 



2 



If X is a real number, then the absolute value of x is denoted \x\ and defined 
to be equal to x when a; > and to —x when x < 0. One can check that 

(1.4) \x + y\ < \x\ + \y\ 

for all real numbers x, y, which is the triangle inequality for the absolute value, 
and that 

(1-5) \xy\^\x\\y\, 

which is to say that the absolute value of a product is equal to the product of 
the individual absolute values. As a consequence of the triangle inequality, one 
can check that 



(1.6) 



\y\ 



<\x-y\ 



for all x,y G H. 

The complex conjugate of a complex number c = a + ib, a, b gR, is denoted 
c and defined by 

(1.7) c = a-ib. 

One can check that 



(1.8) a + p = a + f3 
and 

(1.9) ^ = S/3 

for all a, /? e C. lfc(zC,c — a + ib, with a, 5 G R, then a and b are called the 
real and imaginary parts of c and arc denoted Re c, Im c, respectively, and we 
have that 

(1.10) Rec=^, Imc=^. 

If c = a + i 6 is a complex number, a,b G R, then the absolute value or 
modulus is denoted |c| and defined by 

(1.11) |c| = Va^ + b^. 
Notice that 

(1.12) I Reel, |Imc| < |c| 
and 

(1.13) |cp=cc. 

It follows that the modulus of a product is equal to the product of the moduli, 
so that 

(1.14) |a/?| = |a||/3| 
for a, /3 G C. 

A basic result states that the triangle inequality holds for complex numbers, 

i.e., 

(1.15) |a + /?| < |a| + |/3| 
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for all complex numbers a, (3. Indeed, 

(1.16) |Q! + /3|^ = \a\'^ + 2'R&a'P+\(3\^ 

< |a|2 + 2|a||/3| + |/?|2 

As in the case of real numbers, it follows from the triangle inequality that 

(1.17) ||a|-|/3|| < |a-/3| 
for a, /3 e C. 

Let n be a positive integer. Wc write R" for the space of n-tuples of real 
numbers, also known as n-dimensional Euclidean space. If a; = {x\, . . . ,Xn) 
and y = {y\, . . . ,yn) are elements of R" and f is a real number, then x + y 
and t X arc defined as elements of R" using coordinatewise addition and scalar 
multiplication, so that 

(1.18) x + y= {xi + yi,...,Xn + yn) 

and 

(1.19) tX = {tXi, . . . ,tXn)- 

The inner product of x, y G R" is denoted {x, y) and defined by 

n 

(1.20) {x,y) = Y,xjyj- 
Notice that 

(1.21) {x,y) = {y,x) 

for all x,y G R". If x, x', y e R" and t, t' G R, then 

(1.22) {tx + t' x\y)=t{x,y)+t' {x',y), 

and of course there is a similar linearity property in y by symmetry. 
The Euclidean norm of a; G R" is defined by 



(1.23) I^I=(E^0 



1/2 



This is the same as 

(1.24) |a;|^ = (a;,a;). 

Note that if we identify a complex number c with the element of R^ given by 
(Rec, Imc), then the modulus of c as a complex number is the same as the 
Euclidean norm of the corresponding element of R^. 

An important feature of the inner product on R" is the Cauchy-Schwarz 
inequality, which states that 

(1-25) \{x,y)\<\x\\y\ 
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for x,y & R". To see this, one can observe that for each real number t, 

(1.26) < {x + ty,x + ty) 

= Ixf + 2t{x,y)+t''\yf. 

The Cauchy-Schwarz inequaUty follows by choosing t appropriately. 
For each x £ R" and each t e R we have that 

(1.27) \tx\ = \t\\x\, 

as one can easily check from the definitions. Note that refers to the absolute 
value of the real number t, while \tx\, \x\ refer to the norms of the vectors tx, 
X in R". We also again have the triangle inequality 

(1.28) \x + y\<\x\ + \y\ 

for x,y e R", since 

(1.29) \x + y\'^ = {x + y,x + y) 

= \x\^ + 2{x,y) + \y\^ 

< \x\'+2\x\\y\ + \y\' 

= i\x\ + \y\r- 

Two other important features of the Euclidean norm on R" are the polar- 
ization identity and the parallelogram law. The polarization identity expresses 
the inner product in terms of the norm, through the formula 

(1.30) (^x,y) = ^(^\x + y\'-\x\'-\y\' 

iov x,y G R". The parallelogram law states that 

(1.31) \x + yf + \x-yf = 2\xf + 2\yf, 

which relates the lengths of the sides of a parallelogram to the lengths of the 
diagonals. 



1.2 Linear algebra in 

Let n be a positive integer. In this subsection we review matters concerning 
linear independence, spanning vectors, linear subspaces, and linear mappings in 
R". We begin with a preliminary fact about systems of linear equations. 

Suppose that I and m are positive integers, and that for each j — 1,. . . ,1 
and k = l,...,m we have a real number aj^k- This leads to a system of I 
homogeneous linear equations in m unknowns, namely the equations 

m 

(1.32) ^aj-ferfe=0 

fe=i 
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for j = 1,. .. ,1. Here ri, . . . , are variables which run through the real num- 
bers. 

A basic fact is that if m > I, then this system of equations has a nontrivial 
solution. In other words, in this case there are real numbers ri, . . . ,r„i which 
satisfy the equations above and which are not all equal to 0. This is not difficult 
to show, by using the equations to solve for one variable in terms of the other 
variables and substituting this into the remaining equations to systematically 
reduce the number of equations. 

A collection of vectors vi,. . . ,Vm in is said to be linearly independent if 
for any choices of real numbers n , . . . , we have that 

m 

(1.33) J2'-JVj=0 
in R" implies that 

(1.34) ri = ---rm=0. 

This is equivalent to saying that a vector x in R" can be expressed in at most 
one way as a sum of the form 

m 

(1.35) Y^^^^i 

for some real numbers Si, . . . , s^. Prom the result mentioned in the previous 
paragraph it follows that m is necessarily less than or equal to n for this to 
happen. 

Two vectors v, w in R" are said to be orthogonal if 

(1.36) {v,w)=0, 

and in this case we write v -L w. A collection of vectors vi, {vm is said to be 
orthogonal if -L Vk fov 1 < j k < m, j ^ k. li vi, . . . ,Vm are orthogonal 
vectors in R" such that \vj\ = 1 for each j, then vi, . . . ,Vm is said to be an 
orthonormal collection of vectors in R". 

A collection of nonzero orthogonal vectors vi, . . . , Vm in R" is automatically 
linearly independent. It is enough to check this when wi, . . . ,Vm are orthonor- 
mal, since one can reduce to this case by multiplying the Vj^s by nonzero real 
numbers so that they have norm equal to 1. In this case, if 

m 

(1.37) x = YsjVj 

for some x G R" and real numbers si, . . . , Sm, then 

(1.38) Sj = {x,Vj) 

for each j, so that the coefficients Sj are determined by x. 
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If wi, . . . , is any collection of vectors in R", then a linear combination of 
Wi, . . . ,wi is any vector in R" of the form 



where ri, . . . ,r; are real numbers. Thus a finite collection of vectors in R" is 
linearly independent if each element of R" which can be written as a linear com- 
bination of vectors in the collection can be written as such a linear combination 
in only one way. A collection of vectors Vi,. . . ,Vm in R"* is said to be linearly 
dependent if it is not linearly independent, and this is equivalent to saying that 
one of the Uj's can be written as a linear combination of the other vectors in 
the collection. 

By a linear subspace of R" we mean a subset L of R" such that € L and 
x,y £ L implies that r x + sy G L for all r,s £ R. U wi, ... ,wi is a, collection 
of vectors in R", then the span of wi, . . . , is denoted 

(1.40) span{wi,...,wi) 

and is the subset of R" consisting of all linear combinations of Wi , . . . , w/ . It is 
easy to see that this is always a linear subspace of R". 

As a basic example, define ei, . . . , e„ G R" by saying that ej has jth coor- 
dinate equal to 1 and all other coordinates equal to 0. It is easy to see that 
ei, . . . , e„ is an orthonormal collection of vectors, and is linearly independent in 
particular. Also, the span of ei, . . . , e„ is equal to R". 

If L is a linear subspace of R" and vi,. . . ,Vm are linearly independent vectors 
in R" whose span is equal to L, then we say that vi, . . . ,Vm form a basis for 
L. Thus every linear combination of Vi, . . . ,Vm hes in L in this case, and every 
element of L can be expressed as a linear combination of the vj 's in exactly one 
way. If the vj 's are also orthogonal, or orthonormal, then we say that they form 
an orthogonal basis, or orthonormal basis, of L. 

For example, the vectors ei, . . . , e„ form an orthonormal basis of R". This 
basis is called the standard basis for R". Of course this basis has n elements. 

In general, if i is a linear subspace of R" and vi , . . . , is a basis for L, then 
the dimension of L is defined to be m. One can check that the dimension of a 
subspace is independent of the choice of basis. Indeed, if lOi, . . . , lOp, Vi,. . . ,Vm 
are vectors in R" such that each Vj lies in the span of the tOfe's and the v/s are 
linearly independent, then m < p. 

Let us note that every linear subspace L of R" has a basis. In the case 
where L = {0}, we interpret this as using the "empty basis", and the dimension 
is equal to 0. In general one can get a basis by first choosing a nonzero vector 
in L, and systematically adding vectors to the collection which are not in the 
span of the vectors already selected, until L is equal to the span of the vectors 
selected and one gets a basis. 

If L is a linear subspace of R" and vi, . . . ,Vk is a collection of linearly 
independent vectors in L, then one can add vectors to this collection if necessary 



(1.39) 




fe=i 
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to get a basis for L. This follows from the same kind of argument as the one for 
showing that there is a basis for each linear subspace of R". As a consequence, 
if Li, L2 are linear subspaces of R" such that Li C L2, then the dimension of 
Li is less than or equal to the dimension of L2. 

Now suppose that L is a linear subspace of R" which is equal to the span of 
a collection of vectors wi,. . . , Wm- In this case the dimension of L is less than 
or equal to m. Indeed, either wi, . . . , are linearly independent, and hence a 
basis of L, or one can remove vectors from this collection without changing the 
span to get a subcoUection which is a basis for L. 

Suppose that Li, L2 are linear subspaces of R" such that 

(1.41) LinL2 = {0}. 

Given a collection of linearly independent vectors in Li and another collection of 
linearly independent vectors in ^2- one c:an combine them to get a collection of 
vectors in R" which is also linearly independent. Thus the sum of the dimensions 
of Li and L2 is less than or equal to n. 

For any two nonempty subsets Ai, A2 of R", put 

(1.42) Ai + A2 = {x eK" : X = ai + a2 for some ai G vli, 02 e A2}. 

If Ai, A2 are linear subspaces of R", then so is the sum Ai+ A2. If Ai, A2 are 
linear subspaces of R" such that Ai + A2 = R", then one can check that the 
sum of the dimensions of Ai, A2 is greater than or equal to n. 

By a linear mapping on R" we mean a mapping T from R" to itself such 
that 

(1.43) T{x + y) = T{x)+T{y) 
for all x,y £ R", and 

(1.44) T{rx)=rT{x) 
for all r G R and x G R". In particular, 

(1.45) r(0) = 0. 

The identity mapping, which takes x to itself for each x G R", is denoted /. 

If Ti, T2 are linear mappings on R", and ri, r2 are real numbers, then the 
linear combination ri Ti + r2T'2 is defined as a linear mapping on R", which 
sends x G R" to ri ri(a;)+r2 r2(a;). The composition of Ti , T2 is denoted T10T2 
and defined by 

(1.46) (TioT2)(a;) =Ti(T2(x)) 

for X G R", and is also a linear mapping on R". Of course if T is a linear 
mapping on R", then 

(1.47) ToI^IoT^T. 

Associated to a linear mapping T on R" we have a n x n matrix (tj,fe) of 
real numbers, given by 

(1-48) i,-,fc = (r(efc),e,). 
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where ei, . . . , e„ is the standard basis for R"". This is equivalent to saying that 

n 

(1.49) T(efc)=^t,-fee,-. 

Conversely, if is an n x n matrix of real numbers, then there is a unique 

linear mapping T on R" associated to this matrix in this manner. 

The matrix associated to the identity operator / is the Kronecker delta 
matrix 

(1.50) 5j^k = 1 when j = k 

= when j ^ k. 

If Ti, T2 are linear mappings on R" and {tji.) and {tji.) are the nxn matrices 
associated to Ti, T2, respectively, then the matrix associated to ri Ti + T2 is 
given by (ri tjj, + r2 t'j,k)- The nxn matrix {tj^^) associated to the composition 
T3 = Ti o T2 is given by 

which is the usual notion of "matrix multiplication" . 

If L is any linear subspace of R", then the orthogonal complement of L in 
R" is denoted L-^ and defined by 

(1.52) = {w e R" : w -L ?; for all v G L}. 

It is easy to see that L-^ is also a linear subspace of R". Furthermore, 

(1.53) LnL^ = {0}. 

Suppose that L is a linear subspace of R" and that wi , . . . , is an orthonor- 
mal basis for L. Define a linear mapping Pj, from R" to itself by 

m 

(1.54) Pi(x)=^(x, 

This is called the orthogonal projection of R" onto L. 
One can check that 

(1.55) Pl{x) gL and x- Pl{x) e 

for each x € R". Also, if y\ and 2/2 are two vectors in L such that x — y\,x — 
2/2 G then j/i = 7/2- In other words, Pl{x) is uniquely determined by the 
properties mentioned above, and in particular Pl{x) does not depend on the 
choice of the orthonormal basis v\, . . . ,Vm ^or L. 

If lui , . . . , iwfe is an orthonormal collection of vectors in R" and u is another 
vector in R" which does not lie in span('u;i, . . . , w^), then one can use the 
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orthogonal projection onto span(wi, . . . , Wk) as just described to define a vector 
Wk+i in R" such that 

(1.56) Wi,...,Wk,'Wk+i 
is an orthonormal collection of vectors in R" and 

(1.57) spa,n{wi,...,Wk,u) = spa.n{wi, . . . ,Wk,Wk+i). 

Using this one can show that every linear subspace of R" has an orthonor- 
mal basis. As a result, for each linear subspace of R" there is an orthogonal 

projection of R" onto that subspace. 
If L is a linear subspace of R" , then 

(1.58) LnL^ = {0} and L + L-^ = R", 

and in particular the sum of the dimensions of L and L-^ is equal to n. Let us 
also observe that 

(1.59) {L^)^ = L. 

Indeed, L C (i"*")"*" follows directly from the definition of the orthogonal com- 
plement, and conversely if x G (L^)^, then x — Pl{x) lies in both L"*" and 
(L-"-)-"-, and hence is equal to 0, so that x = Pl{x) G L. 

By a bilinear form we mean a function B{x,y) on R" x R" with values in 
the real numbers which is linear in each of x and y, i.e., 

(1.60) B{r x + r' x', y) = r B{x, y) + r' B{x', y) 
for all r, r' G R and all x, x' , y G R", and 

(1.61) B{x, s y + s' y') = s B{x, y) + s' B{x, y') 

for all s, s' e R and all x, y, y' £ R". If T is a linear mapping on R", then 

(1.62) Bix,y) = {Tix),y) 

defines a bilinear form on R". Conversely, one can check that every bilinear 
form on R" arises in this manner. 

If B(x,y) is a bilinear form on R", then we can get a new bilinear form 
B*{x,y) simply by interchanging the variables, i.e., 

(1.63) B*{x,y) = B{y,x). 

If T is a linear transformation on R", then there is a unique linear transforma- 
tion T* on R", called the adjoint of T, such that 

(1.64) {T*{x),y)^{x,T{y)). 

If (ij,fc), {t*j^k) matrices associated to T, T* , respectively, then 

(1-65) = tkj 
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for j,k = 1, . . . , n. 

If ri, r2 are real numbers and Ti, T2 are linear mappings on R", then 



(1.66) 



{nT2 + r2T2)* =nT^ + r2T^. 



Also, 

(1.67) (Ti 0T2)* = oTi*. 

For any linear transformation T on R", 



(1.68) 



A bilinear form B{x, y) on R" is said to be symmetric if B* = B, a linear 
mapping T on R" is said to be symmetric if T* = T, and an n x n matrix 
(tj^k) is said to be symmetric if tj^k = tk,j- These conditions are related in 
the obvious way, so that a linear transformation is symmetric if and only if 
the corresponding bilinear form is symmetric, and this holds if and only if the 
associated matrix is symmetric. An n x n matrix {tj^k) is said to be diagonal if 
tj,k = when j 7^ fc, which clearly implies that the matrix is symmetric. 

It is easy to see that the adjoint of the identity mapping is equal to itself. If 
Z/ is a linear subspace of R", then the orthogonal projection Pj, of R" onto L 
is symmetric. For if a;, y G R", then 



Let T be a linear mapping on R". The kernel of T is defined to be the set 
of a; e R" such that ^(a;) = 0, and the image is defined to be the set of y e R" 
such that y = T{x) for some x € R". It is easy to see that the kernel and image 
of T are linear subspaces of R". 

It is a simple exercise to show that a linear mapping T on R" is one-to-one, 
which means that T(x') = T{x) implies x' = x for x, x' € R", if and only if the 
kernel of T is trivial, in the sense that it is equal to {0}. One can always restrict 
a linear mapping T to the orthogonal complement of its kernel, and on this 
subspace of R" T is onc-to-onc. Using this one can check that ii vi, . . . ,Vm is a 
basis for the orthogonal complement of the kernel of T, then T{vi), . . . ,T{vm) 
is a basis for the image of T. 

As a result, the dimension of the orthogonal complement of the kernel of T 
is equal to the dimension of the image of T. Hence the sum of the dimensions 
of the kernel of T and the image of T is equal to n. This implies in turn that T 
is one-to-one if and only if T maps R" onto R", which is to say that the image 
of T is aU of R". 

Let us note as well that if y G R" lies in the kernel of the adjoint T* of 
a linear mapping T on R", so that T*{y) = 0, then y lies in the orthogonal 
complement of the image of T. Conversely, if y G R" is in the orthogonal 
complement of the image of T, so that 



(1.69) 



(Pi(x),y) = (Pi(x),Pi(y)) = (x,Pi(y)). 



(1.70) 



(y,T(x))=0 
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for all X e R", then 

(1.71) {T*{v),x)=0 

for all X e R", which implies that T*{y) = 0. In short, the kernel of T* is equal 
to the orthogonal complement of the image of T, and thus the dimensions of 
the kernels of T and T* arc equal to each other. 

If T is a linear transformation on R" and L is a linear subspace of R" , then 
L is said to be invariant under T if T{L) C L, which is to say that T{x) e L 
for all X € L. One can check that if L is invariant under T, then the orthogonal 
complement L-^ of L is invariant under the adjoint T* of T. As a special case, 
if T is a symmetric linear transformation on R" and L is a linear subspace of 
R" which is invariant under L, then L-*- is also invariant under T. 

A linear transformation T on R" is said to be invertible if there is another 
linear transformation S on R" such that 

(1.72) SoT = ToS = I. 

In this case S is called the inverse and is denoted T^^. It is easy to see that 
the inverse is unique when it exists, and that if Ti, T2 are invertible linear 
transformations on R", then Ti o T2 is also invertible, with 

(1.73) (Ti o Ts)-^ = T^^ o Tf 1. 

If T is a linear mapping on R" which is invertible simply as a mapping, which 
is to say that there is a mapping S from R" to itself such that SoT = ToS = I, 
then it is easy to see that S is linear, so that T is invertible as a linear mapping. 
Of course T is invertible as a mapping on R" if it is one-to-one and maps R" 
onto itself. For linear mapping these two conditions are equivalent to each other, 
as noted before. 

A linear transformation T on R" is said to be an orthogonal transformation 

if 

(1.74) {T{x),T{y)) = {x,y) 

for all x,y € R". Because of the polarization identity, this is equivalent to 
saying that 

(1.75) |T(x)| = \x\ 

for all X G R". One can check that this is also equivalent to saying that T is 
invertible, and that 

(1.76) rp-l^rp* 

If T is a linear mapping on R", then the norm of T is denoted ||T|| and 
defined by 

(1.77) ||T|| = sup{|T(a;)| : a; e R", |a;| = 1}. 

It is easy to check that the nonnegative real numbers [^(a;)!, x € R", \x\ = 1, 
are bounded from above, so that the supremum in this definition makes sense. 
To put it another way, the norm ||r|l of T is the nonnegative real number such 
that 

(1.78) |r(x)| < IITII |x| 
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for all X G R", and which is as small as possible. 

Clearly ||T|| = if and only if T = 0, and the norm of the identity operator 
is equal to 1. One can check that 

(1.79) ||rT|| = |r||!r|| 

for every real number r and every linear operator T on R", and that 

(1-80) ||ri + T2||<||ri|| + ||r2|| 

for all linear mappings Ti, T2 on R". Moreover, 
(1-81) ||TioT2||< IIT1IIIIT2II 

for all linear mappings T\, T2 on R". 

An equivalent definition of the norm of a linear transformation T on R" is 

(1.82) ||T|| = sup{|(T(a;),y)| : x,y e R", \x\ = \y\ = 1}. 

Indeed, the right side is less than or equal to ||T|| because 

(1-83) |(T(x),y)|<||r|| 

when x,y & R", = \y\ = 1, by the Cauchy-Schwarz inequality. To get the 
reverse inequality, one can consider the case where y is a, multiple of T{x). 
From this equivalent definition of the norm it follows that 

(1.84) ||T*|| = ||T|| 

for any linear mapping on R". Another basic fact is the "C* identity" 

(1.85) ||T* oTll = ||r||2. 
This can be checked using the simple formula 

(1.86) {{T*oT){x),y) = {T{x),Tiy)). 

If L is a linear subspace of R" which contains a nonzero element, and Pl is 
the orthogonal projection of R" onto L, then 

(1.87) IIP^II = 1. 
Indeed, if x, y are vectors in R" such that x -Ly, then 

(1.88) |x + 2/p = |xp + |j/|2, 
and thus 

(1.89) \Pl{u)\^ < \PLiu)f + \u- Pl{u)\^ = \u\^ 

for all u G R". This implies that \\Pl\\ < 1, and the reverse inequalities hold 
because Pl{u) = u when u £ L. 
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If T is an orthogonal linear transformation on R", then of course 

(1.90) IITII = 1. 

More precisely, an invertible linear mapping T on R" is an orthogonal transfor- 
mation if and only if 

(1.91) IITII = ||T-i|l = l. 

This is easy to check from the definitions. 

A symmetric linear transformation T on R" is said to be nonnegative if 

(1.92) {T{x),x)>0 

for all X e R". More generally, if A, B are symmetric linear transformations on 
R", then we say that A is less than or equal to B, written 

(1.93) A<B, 

if i? — ^ is nonnegative. For any symmetric linear transformation T on R", we 
have that 

(1.94) - ||T||7<T< ||T||7. 

If T is a nonnegative symmetric linear operator on R", then 

(1.95) mx),y)\ < {T{x),x)''^ {T{y),y)"^. 

This can be checked in the same manner as for the Cauchy-Schwarz inequality. 

In particular, T{x) — when (T(x), x) = 0. 

If T is a symmetric linear transformation on R", then 

(1.96) (r(a;),a;)>0 

for all X G R" with .x if and only if T is invertible. Indeed, if this condition 
holds, then the kernel of T contains only the zero vector, and T is invertible. 
Conversely, if T is invertible, then T{x) ^ for all x e R", and hence the 
inner product of T{x) and x is nonzero by the result mentioned in the previous 
paragraph. 

1.3 Sequences and series 

Fix a positive; integer n, and let {xj}'jli be a sequence of points in R". This 
sequence is said to converge to a point x G R" if for every e > there is an 
integer N = N{e) so that 

(1.97) |xj -x\<e for all j > N. 

In this case x is said to be the limit of the sequence {xj}JLi, which one can 
easily show is unique, and one writes 

(1.98) lim Xj = X. 

j—>oo 



14 



Sometimes it will be convenient to consider a sequence {a;j}°^(j with a dif- 
ferent starting point a £ Z, and the same basic notions and results apply just as 
well. The same definition also applies to sequences of complex numbers, using 
the standard identification of C with R^. In this regard let us note that a se- 
quence {xj }j in R" converges to a point x G R" if and only if for each integer I 
such that 1 < I <n the sequence of Zth components of the Xj 's converges to the 
Ith component of x as a seqiience of real mimbcrs, and in particular a sequence 
{zj}j of complex numbers converges to a complex number z if and only if the 
real and imaginary parts of the Zj 's converge to the real and imaginary parts of 
z, respectively. 

Suppose that {xj}^!^, {yj}j^i are sequences in R" which converge to x,y £ 
R", respectively. A standard result which is easy to verify states that {xj + 
yj}j^i converges to x + y. In particular this applies to sequences of real and 
complex numbers. 

Now suppose that {zj , {wj are sequences of complex numbers which 
converge to the complex numbers z, w, respectively. One can then show that 
the sequence of products {zj Wj}'j^i converges to the product zw of the limits of 
the individual sequences. This includes products of sequences of real numbers 
as a special case, and for a pair of sequences {xj}"^-^^, {yj}°^i in R" which 
converge to x,y € R", respectively, one has that the sequence of inner products 
{{xj,yj)}'j^i converges to the inner product {x,y) of the limits. 

One way to deal with sequences of products is to consider separately the 
cases where one of the sequences is constant, and where one of the sequences 
converges to 0. If one of the sequences tends to 0, then it is sufficient for the 
other sequence to be bounded in order for the sequence of products to also 
tend to 0. Let us note as well that if {zj}jli is a sequence of nonzero complex 
numbers which converges to a nonzero complex number z, then the sequence of 
reciprocals {^/zj}"^-^ converges to the reciprocal of the limit, 1/z. 

Sequences of linear transformations on R" can be converted into sequences 
of n X n matrices, which can then be reinterpreted as sequences in R™ with 
m = n^, so that convergence can be defined as before. If {Sj}°^i, {Tj}°^i 
are sequences of linear transformations on R" which converge to the linear 
transformations S, T, respectively, then {Sj+Tj}'^-^ and {SjoT,}^^ converge 
to S* + T and S o T, respectively. If {Tj}°^^ is a sequence of invertible linear 
transformations on R" which converges to the invertible linear transformation 
T, then {T^^}jZi converges to T~^, as one can show using "Cramer's rule" for 
expressing the inverse of a matrix in terms of determinants. 

Suppose that {xj}j^i is a sequence of real numbers. We say that this se- 
quence is monotone increasing if xj+i > xj for all j. Similarly, {xj}^!^ is said 
to be monotone decreasing if .Xj+i < Xj for all j. 

Let {xj}°°^i be a monotone increasing sequence of real numbers which is 
bounded from above, which is to say that there is a real number b such that 
Xj < b for all j. In this case {xj}jli converges to a real number x, with 

(1.99) X = sup{xj : i G Z+}. 

In the same way, if {y^jj^i is a monotone decreasing sequence of real numbers 
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which is bounded from below, then {yj}'j^i converges to y = infj>i yj. 

A sequence {xjjJLi in R" is said to be a Cauchy sequence if for every e > 
there is a positive integer N such that 

(1.100) \xj -Xk\<e for aU j, k>N. 

It is easy to see that a convergent sequence is a Cauchy sequence, and in fact the 
converse holds too, which is to say that every Cauchy sequence converges. To 
show this it is enough to consider a Cauchy sequence {xj}°^i of real mmibers, 
and the limit of this sequence is the same as the limit of the monotone increasing 
sequence {x'j}°^^ defined by Xj = infj>j Xj. 

Now let J2jLo % be a series of complex numbers. Associated to this series 
is the sequence of partial sums 

n 

(1.101) ^a,-, 

where n runs through the nonnegative integers. If the sequence of partial sums 
converges to a complex number A, then we say that the series YlJLo converges 
to A. 

Suppose that J2jLo '^j ' S jlo series of complex numbers which 

converge to the complex numbers A, B, respectively. If a, (3 are complex num- 
bers, then the series 

oo 

(1.102) "Y^aaj+phj 

converges to aA + ^B. This follows from the corresponding statement for 
sequences, since the partial sums of the new series arc equal to the same kind 
of linear combinations of the partial sums of the original series. 

The Cauchy criterion for convergence of a series states that a series X] jlo 
of complex numbers converges if and only if for every e > there is a positive 
integer N such that 



(1.103) 



< e 



for all TO, n such that n> m> N. Indeed, this condition is equivalent to saying 

that the sequence of partial sums of the series X]^o '^3 ^ Cauchy sequence. A 
consequence of the Cauchy criterion is the comparison test for convergence of a 
series, which says that if Y^JLo a series of complex numbers and Y^JLo 
is a scries of nonnegative real numbers which converges such that \aj\ < bj for 
all j, then X)^o % converges. 

Suppose that X^^lo ^ series whose terms are nonnegative real numbers. 
Then the scries converges if and only if the sequence of partial simis is bounded 
from above. This is because the sequence of partial sums is a monotone increas- 
ing sequence of real numbers. 

The results mentioned in the preceeding two paragraphs are quite simple, 
and also quite useful. In short, the question of convergence of a series can often 
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be reduced to finding an upper bound for some collection of nonnegative real 
numbers. This is often an easier task. 

Let us note a basic necessary condition for a series to converge. Namely, if 
Ylj=n '^j ^ scries of complex numbers which converges, then {o jlj^o converges 
to 0. This is easy to derive from the Cauchy criterion for convergence. 

A series X^^o ^" converges absolutely if X^^Iq converges, which implies 
convergence of the original scries in the ordinary sense. Absolute convergence is 
a more stable kind of convergence. For example, if a series converges absolutely, 
then one can rearrange the terms and still get a series which converges, and 
with the same s\im. 

A famous test for convergence of series is called the Cauchy Condensation 
Test, and it says the following. If Y^^=i ^ series of nonnegative real numbers 
such that the sequence of bj's is monotone decreasing, then J2'jLi converges 
if and only if the "condensed" series Yl'k=o '^^^2'' converges. To prove this, one 
shows that the partial sums of each can be bounded in terms of the other. 

Before describing some examples, let us review some auxiliary facts. First, 
the binomial theorem states that for each positive integer n, 

(1.104) {x + yY = Y,(^\x^y^-\ 



where 
(1.105) 



m! = 1 • 2 • • • m, 0! = 1. This can be proved by induction, with the n = 1 case 
being trivial, and let us note also that the binomial coefficients (^) are always 
integers. 

Now let o be a positive real number, and observe that 

(1.106) (l + a)">l + na 

for all positive integers n, by the binomial theorem. Thus 

(1.107) lim (1 + a)-" = 0, 

and hence 

(1.108) lim z"- = 

n — ^oo 

when z is a complex number with \z\ < 1. More generally, if z is a complex 
number with < 1 and fc is a positive integer, then 



(1.109) lim n'^z" = 0. 

n — ^oo 



If z is a complex number with \z\ < 1, then the geometric series Y1'^=^qZ" 
converges, and 

°° 1 

(1.110) E^' = T3T- 

3=0 
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Indeed, 

n 

(1.111) (1 - ^) = 1 - 

3=0 

for each nonnegative integer n, so that 



(1.112) lim V 



hm 



1 - z 



n+l 



_ n^oo 1 — Z 1 — Z 

Of course the scries docs not converge when \z\ > 1, since \z\" > 1 then for all 
n, and the terms of the series do not tend to 0. 

Now suppose that p is a positive real number, and consider the series 



oo -J 



, nP 

n—l 

One can show that this series converges if and only if p > 1, using the Cauchy 
Condensation Test. More precisely, this reduces the question to one for geomet- 
ric scries. 

However, the series 

- (-1)" 



(1.114) y 

converges for all positive real numbers. This can be viewed as a special case of 
Leibniz' alternating series test, which states that a series of the form 

oo 

(1.115) 

converges if {6„}^j is a monotone decreasing sequence of real numbers which 

converges to 0. One can verify this test using the Cauchy criterion, by estimating 
the absolute value of sums of the form J2n=pi~^)" ^" terms of the absolute 
value of the first term bp. 

More generally, a series of the form 

oo 

(1.116) E'^"^" 

converges if the sequence of partial sums Ap = YX,=i of the a„'s is bounded 

and if {&ri}J^]^ is a monotone decreasing sequence of real numbers which con- 
verges to 0. Indeed, if we set Aq = for convenience, then for each positive 
integer p, 

p p 
(1-117) E 

= J^i^ri- An-l)bn 
n=l n=l 
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p 

^ ' A-n bn — 
n=l 


P 
n=l 


P 


P-1 


^ ^ j4„ bn — 




ri = l 




p-1 




J2 ^» ■ 


-bn+i) + Apbp 



One can check that J2'^=i {bn — ^n+i) converges absolutely, since 

oo 

(1-118) Y.^bn-bn+l) 

n=l 

is a convergent "telescoping series" with nonnegative terms, and hence that 

«n bn converges. 

For example, suppose that 2; is a complex number such that |2;| = 1 and 
z ^ 1, and consider a„ = 2;". In this case we have that 

P r _ rP 

(1.119) E-" = ^' 

n=l 

which is a bounded sequence of complex numbers as p runs through the positive 
integers. Thus we obtain that 

00 

(1.120) 

n=l 

converges when z is a complex number such that \z\ = 1 and z ^ 1 and {bn}^—i 
is a monotone decreasing sequence of real numbers which converges to 0. 

1.4 QuEirternions and Clifford algebras 

The quarternions are denoted H, and a quarternion can be written as a + bi + 
cj + dk, where a, b, c, d are real numbers, and i, j, k satisfy the conditions 

(1.121) f=f = k'^ = -l 
and 

(1.122) k = ij = -ji. 

To be more precise, the quarternions arc an algebra over the real numbers, with 
coordinatewise addition, the real number 1 being the multiplicative identity el- 
ement for the whole algebra, and real numbers commuting multiplicatively with 
other quarternions even if multiplication in H is not commutative in general. 
Multiplication in H does satisfy the usual associative and distributive rules, and 

(1.123) ik = —ki = —k, jk = kj = —i 
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in particular. 

Suppose that x = a + bi + cj+dk is a quarternion. The conjugate of x is 
denoted x* and defined by 

(1.124) a* =a-bi~ck-dk. 
li X, y are quarternions, then {x*)* = x, 

(1.125) {x + yy=x*+y*, 
and 

(1.126) {xyr^y*x*. 

The modulus of a quarternion x = a + bi + cj + dk is denoted |a;| and defined 

by 

(1.127) \x\ = {a^ + b^ + c'' + d^)^/\ 
One can check that 

(1.128) \x\^ ^ XX* = x* X, 

and in particular it follows that a nonzero quarternion x is invertible, with 
= \x\~^ x*. Moreover, 

(1.129) \xy\ = \x\\y\ 

for x, 2/ G H. 

If x = a + bi + cj + dk is a quarternion, with a, b,c,d€ R, as usual, then a 
is called the real part of x. Of course the real part of x is equal to the real part 
of X* . One can check that if x, y are quarternions, then the real part of a; y is 
equal to the real part of yx. 

We can identify a quarternion x = a + bi + cj + dk, a,b,c,d gH, with the 
element (a, 5, c, d) of R'*. This identification respects the operations of addition 
and "scalar" multiplication by real numbers on H and R'*, and the modulus of 
a quarternion corresponds exactly to the standard Euclidean norm on R*. If x, 
y are quarternions, then the real part of x y* is the same as the standard inner 
product of the cooresponding elements of R^. 

Suppose that a, (3 are quarternions such that \a\ = \(3\ = 1, and consider 
the mapping 

(1.130) x^axf] 

from H onto itself. If we identify H with R** as in the previous paragraph, 
then this cooresponds to an orthogonal linear transformation on R^. This is 
analogous to the fact that if ^ is a complex number such that |^| = 1, then the 

mapping 

(1.131) z^Oz 

from C onto itself corresponds to an orthogonal linear transformation on R^ 
using the standard identification of C with R^. 

A quarternion w is said to be imaginary if its real part is equal to 0, which 
is equivalent to 

(1.132) w* = -w, 
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and to 
(1.133) 



w^ = -\w\\ 



If a is a quarternion and w is an imaginary quarternion, then aw a* is also an 
imaginary quarternion. If a is a quarternion with \a\ = 1, then the mapping 

(1.134) w 1-^ a w a* = a w 

on imaginary quarternions corresponds, under the natural identification of imag- 
inary quarternions with elements of R^, to an orthogonal Hnear transformation 
on R^. 

Now let n be a positive integer, and let C(n) denote the Clifford algebra 
with n generators ei, . . . , e„. To be more precise, C{n) is an algebra over the 
real numbers which contains a copy of the real numbers such that the real 
number 1 the multiplicative identity element for the whole Clifford algebra and 
of course real numbers commute multiplicatively with all other elements of the 
Clifford algebra. The generators ei, . . . , e„ satisfy the relations 

(1.135) e| = -1 
for j = 1, . . . , n and 

(1.136) ejek = -ekej 
when j ^ k. 

As an algebra over the real numbers, C{n) is a vector space over the real 
numbers. As such it has finite dimension. The basic reason for this is that any 
product of gj's can be reduced to one in which each e; appears at most once. 

In fact, the dimension of C(n) as a vector space over the real numbers is 
equal to 2". For let / be a subset of the set {1, . . . , n}, where I consists of 

(1.137) ii<i2<---<ii 

for some I, < I < n, and let ej denote the product 

(1.138) • • • Ci, , 

where this is interpreted as being equal to 1 when 1 = 0. The family of e/'s, 
where / runs through all subsets of {1, . . . ,n}, forms a basis for C{n), which 
is to say that every element of C(n) can be expressed in a unique manner as a 
linear combination of the e/'s. 

One can define C(n) when n = to be just the real numbers themselves. 
When n = 1, the Clifford algebra C(n) is equivalent to the complex numbers. 
When n = 2, the Clifford algebra is equivalent to the quarternions. 

Unlike the real numbers, complex numbers, or quarternions, it is not the 
case in general that every nonzero clement of C{n) has a multiplicative inverse. 
However, suppose that x €C{n) is of the form 

n 

(1.139) X = Xo + y~^a;j ej, 
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where xq, a;i, . . . , a;„ e R. If a; 7^ 0, then x is invertible in C{n), with 



(1.140) 



,T|| 



^0 + Sj=i 



2 Harmonic and holomorphic functions, 1 

If X £ R" and r is a positive real number, then we define B(x, r), the open ball 
with center x and radius r in R" by 

(2.1) B(x,r)={2/eR":|y-a;|<r}. 

The closed ball with center x and radius r is denoted B{x,r) and defined by 

(2.2) B{x,r) = {y&W' ■.\y-x\<r}. 

A subset E of R" is said to be hounded if it is contained in a ball. 

A subset U of R" is said to be open if for every x € U there is a positive 
real number r such that 

(2.3) B{x, r) C U. 

The empty set and R" itself are automatically open subsets of R", and one can 
check that open balls in R" are open subsets, using the triangle inequality for 
the Euclidean distance \v — iu\, v,iu € R". The interior of a subset A of R" is 
denoted A° and is defined to be the set of x & A for which there is a positive 
real number r such that B{x,r) C A, and one can check that the interior of A 
is always an open subset of R". 

If is a subset of R" and p is a point in R", then p is said to be a limit 
point of E if for each positive real number r > there is an a; G such that 
X ^ p and 

(2.4) \x - p\ < r. 

A subset E of R" is said to be closed if every point in R" which is a limit point 
of E is also an element of E. This is equivalent to saying that every sequence 
{xj}j^-i of points in E which converges to some point x G R" has its limit x in 
E. 

The empty set and R" itself arc automatically closed subsets of R", and one 
can check that a closed ball in R" is a closed subset of R". If is a subset of 
R", then the closure of E is denoted E and is defined to be the union of E and 
the set of limit points of E, and it is not difficult to show that the closure of a 
subset of R" is always a closed subset of R". If A is a subset of R", then A is 
open if and only if the complement R"\j4 of A in R", consisting of the points 
in R" not in A, is a closed subset of R". 
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2.1 Some differential operators 

Let ?7 be a nonempty open subset of R", and let f{x) be a twice continuously- 
difFerentiable real or complex-valued function on U. Put 

j=l J 

This is called the Laplacian of /. 

Now suppose that n = 2, and identify with C. Writing z = x + iy, 
x,y € R, for an element of C, if f{z) is a continuously-differentiable complex- 
valued function on a nonempty open subset U oi C, then we put 

and 

If f{z) is twice continuously-differentiable on U, then we have that 

(2.8) AA/= _^_^/= Ia/. 

dz dz dzdz 4 

For general n again, assume that f{x) is a continuously-differentiable func- 
tion on a nonempty open subset U of R" with values in the Clifford algebra 
C(n). Define the corresponding left and right Dirac operators acting on / by 

(2.9) ^^/ = E^.a^/ 
and 

(2-10) ^«/ = Ea^/^^' 

where ei, . . . , e„ are the usual multiplicative basis for C{n), i.e., we either multi- 
ply the first derivatives of / by the e^'s on the left or on the right, respectively. 
If / is twice continuously-differentiable, then we have that 

(2.11) Vlf = Vy = -A/, 

where the Laplacian of /, A/, is defined for a Clifford-valued function in the 
same way as before, which amounts to applying the Laplacian for real-valued 
functions to the components of /. 

Let us note that there are some natural variants of this, where for instance 

one considers functions on R"+-'^ with values in C{n), and one uses 1 as well as 
ei, . . . , e„ for defining the first-order differential operators. This type of set-up 
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is more directly analogous to the one for complex- valued functions on C. For 
simplicity we shall use the version described in the previous paragraph. 

If p is a point in R", r is a positive real number, and f{x) is a continuous 
function on the closed ball B{p, r) in R", then we write averageg^^ ^-^ f for the 

average of / over B{p, r), which is the integral of / over B{p, r) divided by the 
volume of B{p^ r). For j, fc = 1, . . . , n, we have that 

(2.12) average^^p^^) {xj - Pjf 
and thus 

(2.13) average5(p_^) {xj - pjf 
Also, 

(2.14) avorage;g(p_^j|a; 

as one can see by reducing to the case where p = and using polar coordinates. 

Let U he B. nonempty open subset of R". Suppose that f{x) is a twice 
continuously-differentiable function on U, and that p is an element of U. One 
can express the Laplacian of / at p as 

n + 2 / \ 

(2.15) A/(p) = lim (^average5(^^,)/ - /(p) j . 

This follows from the Taylor approximation of / at p of degree 2. 

A twice continuously-differentiable function / on a nonempty open subset 
U of R" is said to be harmonic if 

(2.16) A/(p) = 
for all p ^ U. This is equivalent to saying that 

(2.17) lim 1 (^average^(^_,)/ - f{p)j = 

for all p G C/. A stronger "mean value property" of harmonic functions will be 
described in the next subsection. 

A continuously-differentiable complex- valued function f{z) on an open sub- 
set ?7 7^ of C is said to be holomorphic if 

(2.18) |/ = 

at every point in U. The derivative f'{z) of a holomorphic function f{z) on U 
is defined by 

(2.19) fiz) = ^J{z). 

If we identify C with R^ in the usual way, then the difii'erential of a continuously- 
differentiable complex-valued function f{z) on {/ at a point p in {/ is a linear 



= average^(p ,^j(a;fc -pfc) , 
= iaverage5(p_^)|x-p|2. 



|2 2 
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mapping from to itself, and in fact {d/&z)f vanishes at p if and only if the 
differential of / at p is given by multiplication by a complex number on C, and 
that complex number is equal to {d/dz)f at p. 

A continuously-diffcrcntiablc C(n)-valucd function f{x) on a nonempty open 
subset U of R" is said to be left or right Clifford holomorphic, respectively, if 

(2.20) VLfix) = 

or 

(2.21) Vi,Jix) = 

for all X € U, respectively. These two conditions are not equivalent in general, 
because of noncommutativity, although there are functions which are both left 
and right Clifford holomorphic. These conditions do not in general imply that 
the differential of / is given in terms of Clifford multiplication, which turns out 
to be a very restrictive condition. 

If a complex or Clifford holorniorphic function is twice continuously differen- 
tiable, then it is also harmonic, in the sense that its components are harmonic, 
as one can see from the identities relating the Laplacian to the various first-order 
differential opc;rators mentioned before. Wc shall sec in the next subsection that 
the condition of twice continuous-differentiability actually holds automatically 
for complex and Clifford holomorphic functions. Note that when n = 1, a 
continuously-diffcrcntiablc function on an open interval (a, b) is "Clifford holo- 
morphic" if and only if its derivative is equal to on the interval, so that the 
function is constant, and a twice continuously-differentiable function on (a, b) is 
harmonic if and only if its second derivative is equal to on the interval, which 
is to say that the function is of the form ax + f3. 

If h{z) is a twice continuously-differentiable function on a nonempty open 
subset J7 of C which is harmonic, then 

(2.22) / = §-h 

is a continuously-differentiable function on U which is holomorphic. If h{x) is 
a twice continuously-differentiable C(n)-valued function on a nonempty open 
subset U of R" which is harmonic, then the functions 

(2.23) .fL = VLh and fR = VRh 

are continuously-differentiable C(n)-valued functions on U which are left and 
right Clifford holomorphic, respectively. If h{x) happens to be real-valued, then 

2.2 Weak derivatives, 1 

Fix a positive integer n, and let U he a nonempty open subset of R". Also let </> 
be a real, complex, or C(n)-valued function on U. We say that ^ has restricted 
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support in U if the set oi x £U such that (j){x) ^ is bounded, and if there is 
a positive real number rj such that 



(2.24) 



x-y\>r} 



when X eU, y e R", and <j){x) ^ 0. 

The sum of two functions with restricted support in U has rerstricted sup- 
port in U, and the product of two functions on U, where at least one of the 
two functions has restricted support in U, also has restricted support in U. We 
shall mostly be concerned with continuous functions on U which have restricted 
support, and which then have nice properties in general, such as being bounded, 
because of basic results from advanced calculus. Notice that if (f){x) is a con- 
tinuous function on U with restricted support, and if we extend (j){x) to all of 
R" by setting (j){x) = when x £ R"\[/, then this extension is a continuous 
function on R" with restricted support in R". 

Suppose that (f){x) be a continuous function on U with restricted support. 
The usual Riemann integral of ^{x) on U is then defined, which we can denote 



As in the previous paragraph, we can just as well think of <j){x) as a continuous 
function on R" with restricted support, and this integral as an integral over 
R", even if only a bounded region in R" is really involved in the integral. 
Let h{x) be a continuous function on U. We say that h is weakly harmonic 



for all twice continuously-differentiable functions (p onU with restricted support 
in U. Notice that if h is twice continuously-differentiable and harmonic in the 
usual sense, then h is harmonic in the weak sense, by integration by parts. 

Let us consider for a moment the special case of a function 0(x) on R" which 
is of the form where p{t) is a twice continuously-differentiable function 

on the real line. By the chain rule, 



Of course we are interested in functions with restricted support, and so let 
us assume that there is a positive real number r such that p{t) = when t>r. 
For simplicity let us assume also that there is a positive real number e such that 
p{t) is constant for t < e. In particular, p'{t), p"{t), and a{t) are then equal to 
for t < e. 

By integration by parts, we automatically have 



(2.25) 




if 



(2.26) 




(2.27) A<P{x)=4\x\^p"i\x\^)+2np'{\x\^). 
In other words, Acj) is of the form cr(|a;p), where 

(2.28) a{t) =4:tp"{t) + 2np'{t). 



(2.29) 




A(j){x) dx = 
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which is equivalent to 

(2.30) / a{t'^)e-Ut = 0, 

Jo 

and to 

/■oo 

(2.31) / CT(i)i("/2)-i^^ = o. 
Define a function 6{u) on R by 

(2.32) 0{u) = Jw""/^ £ a{t) dt. 

Thus 0{u) = when u < e, and also when u> r. 
By construction, 

(2.33) 4 w6l' (u) + 2n6'(M) = cr(M). 

This is the same equation that p' satisfies. It follows that 0{u) = p'{u) for all 

MGR. 

Conversely, suppose that (7{t) is a continuous function on the real line which 
satisifies the conditions mentioned above, i.e., a{t) = when t < e and when 
t > r, where e, r are positive real numbers, and the same integral as before is 
equal to 0. Then we can define 9 again using the previous formula, so that is 
a continuous function on R such that 0(u) = when u < e and when t > r in 
particular. We can define p on the real line so that p' = and p{u) = when 
u > r, and we also have that p{u) is constant on the set of u e R such that 
u < e. 

In other words, we can reverse the process and start with a continuous 

function a on the real line satisfying the conditions in the previous paragraph, 
and obtain a twice continuously-differentiable function p on the real line such 
that p{u) is constant when u < e and equal to when u>r. Prom the function 
p on the real line we get a function (f){x) = p{\x\'^) on R". We want to use these 
functions as test functions for the weak harmonicity property. 

Actually, we would like to use translates of such a function </>. For each 
p e R", put 

(2.34) cl>p{x) = c^{x-p)=p{\x-p\^). 
Thus 

(2.35) A(t>p{x)=a{\x-p\^). 

Let us digress a moment with a useful definition. If A is a nonempty subset 
of R" and x is a point in R", then the distance from a; to A is denoted dist(a;, ^4) 
and defined by 

(2.36) dist(x, A) = ini{\x - a\ : a G A}. 

One can check that dist(a::, A) = if and only if x lies in the closure of A, and 
in particular this holds if and only if x £ A when A is a closed subset of R". 

Now suppose that p is an element of U, and that r is a positive real number 
which satisfies 

(2.37) < dist(p,R"\l7) 



27 



if ?7 is a proper subset of R", and otherwise is arbitrary if [/ = R". Suppose 
that a{t) is a continuous function on the real Une such that a{t) = when t < e 
for some e > 0, a{t) = when t > r for the choice of r being used now, and 
a(t) satisfies the integral condition discussed before. Thus we get associated 
functions 9, p on the real line and ^p(a;) = p{\x — pp) on R", and (/>j, has 
restricted support in U because of the constraint on r. 

Let us apply this to our weakly harmonic function /i(a;) on U . Namely, we 
have that 

(2.38) / (A(/)p(a;))/i(a;)dx = 0, 

JU 

which is to say that 

(2.39) / a{\x-vf)h{x)dx = ^. 

Ju 

We want to use this with interesting choices of a. 

We can reformulate this condition by saying that if ^{x) is a continuous 
radial function on R", so that il){x) is actually a function of a is a positive 
real number such that ^(2;) = when |a;| > a, p is a point in U such that 

(2.40) dist(p,R"\?7) > a, 
and V' satisfies the integral condition 

(2.41) / tp{x) dx = 0, 
then 

(2.42) / tl}{x-p)h{x)dx = Q. 

Ju 

In other words, '4>{x) corresponds to a{\x\^). We have dropped the condition 
that i){x) be equal to when |a:| is sufficiently small, because one can reduce to 
this case through an approximation argument. 

Let us reformulate this again as follows. Suppose that &i(x), 62(2;) are con- 
tinuous radial functions on R", a is a positive real number such that hi{x) = 
when \x\ > a, i = 1, 2, p is an element of R" such that the distance from p to 
R"\f/ is greater than a, and 

(2.43) / bi{x)dx= [ b2{x)dx = l. 
Jr" Jh" 

Then 

(2.44) / bi{x — p) h{x) dx = / b2{x — p) h{x) dx, 
Ju Ju 

which follows from the statement in the previous paragraph by taking iIj{x) = 
bi{x) - b2{x). 
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This impUcs that if b{x) is a continuous radial function on R", a is a positive 
real number such that b{x) = when \x\ > a,pis an element of U whose distance 
to R"\f7 is greater than a, and 

(2.45) / b{x) dx = 1, 
then 

(2.46) / b{x-p)h{x)dx = h{p). 

Ju 

To get this from the previous statement, one can take bi{x) = b{x), and choose 
b2{x) so that it is concentrated as near to p as one wants. Because h is contin- 
uous, the integral of 62 times h is approximately equal to h{p), and in the limit 
we get the desired formula. 

In fact, if p is a point in U and t is a positive real number such that 

(2.47) t<dist(p,R"\Z7), 
then 

(2.48) averageg(p = h{p). 

This is called the mean value property for the continuous function h on U. It is 

easy to derive this from the previous statement by approximating the average 
of /i on a closed ball in U by integrals against continuous radial functions, 
and conversely one can derive the previous assertion for integrals of h against 
continuous radial functions from this one about averages on closed balls. 

To summarize, a continuous function h{x) on a nonempty open subset of R" 
which is weakly harmonic also satisfies the aforementioned mean value property. 



2.3 Weak derivatives, 2 

Let [/ be a nonempty open subset of R", and let h he a. continuous function 
on U which satisfies the mean value property. Let b{x) be a continuous radial 
function on R" and r a positive real number such that b{x) = when |a;| > r 
and 

(2.49) / b{x)dx = l. 

Jr" 

On the set 

(2.50) Ur = {peU : dist(p, R"\[/) > r}, 
which one can check is an open subset of U, we have that 

(2.51) h{p) = I b{x- p) h{x) dx. 

Ju 

We may as well take b{x) to be as smooth as wc like here, and so we choose 
b{x) so that it is continuously differentiable of all orders. By differentiating 
under the integral sign, it follows that h{x) is in fact continuously differentiable 
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of all orders on Ur- This works for all r > 0, and therefore h{x) is continuously 
differentiable of all orders on U. 

Thus a continuous function which is weakly harmonic is also continuously 
differentiable of all orders, and twice continuously-diffcrentiable in particular. 
Such a function is harmonic in the usual sense. Similarly, a continuous function 
which satisfies the mean value property is harmonic. 

Recall that an open subset U of R" is said to be connected if for every 
pair of points x,y G U there is a continuous path in U that goes from x to 
y; equivalently, U is not connected if it can be expressed as the union of two 
disjoint nonempty open subsets of R". If h is a continuous real- valued function 
on a nonempty connected open subset U of R" which satisfies the mean value 
property, and if p is a point in U such that 

(2.52) h{x) < hip) 
for all X €U, then one can check that 

(2.53) h{x) = hip) 

for all x €U, because {x G U : hix) < hip)} is automatically an open subset of 
U since h is continuous, while {x & U : hix) = hip)} is a nonempty open subset 
of U under the present conditions. For the same reasons, if g is a point in U 
such that 

(2.54) hix) > hiq) 
for all a: e U, then 

(2.55) hix) = hiq) 
for all x gU. 



2.4 Weeik derivatives, 3 

Let U he a, nonempty open subset of C, and let fiz) be a continuous complex- 
valued function on U. We say that fiz) is weakly holomorphic if 

(2.56) / (-lzi>iz)]fiz)dz = 



for all continuously differentiable complex- valued functions tpiz) on U with re- 
stricted support in U. If /(z) is continuously differentiable and holomorphic, 
then fiz) is weakly holomorphic, by integration by parts. 

Suppose that fiz) is weakly holomorphic on U. We can choose the test 
function ipiz) above to be of the form 

(2.57) i^iz) = ^cj^iz), 

oz 

where (piz) is a twice continuously differentiable function with restricted support 
in U, and it follows that 



(2.58) iA4>iz))fiz)dz = 



( iAm)fiz), 
Ju 
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so that f{z) is weakly harmonic. Thus wc obtain that f{z) is continuously 
differentiable of all orders on U, and that f{z) is holomorphic in the usual 
sense. 

Now let [/ be a nonempty open subset of R". and let f{x) be a continuous 
C(n)-valued function on U. We say that f{x) is weakly left Clifford holomorphic 
if 

(2.59) / VRiPi{x)f{x)dx = 

Ju 

for all continuously-differentiable C(n)-valued functions Vi with restricted sup- 
port in U, and we say that f{x) is weakly right Clifford holomorphic if 

(2.60) / f{x)VLMx)dx = 

Ju 

for all continuously-differentiable C(n)-valucd functions tp2 with restricted sup- 
port in U. Again one can choose tpi, ^'2 to be of the form Vj^cj), V^cj) for a twice 
continuously-differentiable function (j) with restricted support in U to obtain 
that a wc;akly left or right Clifford holomorphic function is weakly harmonic, 
and hence harmonic, and thus continuously differentiable of all orders and left 
or right Clifford holomorphic function in the usual sense. 

3 Harmonic and holomorphic functions, 2 

3.1 Poisson kernels, 1 

Fix a positive integer n. By a multi-index we mean an n-tuple 

(3.1) a= (ai,...,a„) 

where each is a nonncgativc integer. In this case the degree of the multi-index 
a is denoted d{a) and defined by 

(3.2) d{a) = ai + --- + an. 

Let a = (qi. . . . , «„) be a multi-index. The corresponding monomial is the 
function of a; = [xi, . . . ,Xn) € R" defined by 

(3.3) 

The degree of the monomial x" is defined to be the degree of the multi-index a. 

By a polynomial on R" we mean a function which is a finite linear combi- 
nation of monomials. A harmonic polynomial is simply a polynomial which is 
harmonic. For a polynomial p{x) on R", the Laplacian Ap(a;) can be expressed 
algeraically, and thus the property of being harmonic can also be expressed 
algebraically. 
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A polynomial p{x) on R" is said to be hom,ogeneous of degree d, where d is a 
nonnegative integer, if p(a;) is a finite linear combination of monomials x" with 
d{a) = d. This is equivalent to saying that 

(3.4) p{tx) = t'^p{x) 

for alH G R and all x G R". Every polynomial on R" can be decomposed into 
homogeneous parts in a simple way. 

If p{x) is a homogeneous polynomial on R" of degree rf, then the Lapacian 
Ap{x) of p{x) is equal to if d = 0,1 and is a homogeneous polynomial of degree 
d—2 when d>2. If h(x) is a harmonic polynomial on R", then its homogeneous 
components are also harmonic. If q(x) is a homogeneous polynomial of degree 
i, then q{x) is a homogeneous polynomial of degree £ + 2. 

Suppose that p{x) is a homogeneous polynomial on R" of degree d. Then 

" d 

(3.5) P(^) = ^P(^)' 

j=i 1 

which is known as "Euler's identity" . This follows by computing the derivative 
oip{tx) in t and setting t = 1, where the derivative is computed using the chain 
rule, and by using homogeneity to write p{t, x) as p{x) and differentiating t''' 
directly. 

If q{x) is a polynomial on R", then 

n Q 

(3.6) \{\x\^ q{x)) = 2nq{x) + 4 ^ x,- — q{x) + \x\'' Aq{x). 

Assuming also that q{x) is homogeneous of degree i, it follows that 

(3.7) A(|a;|2 q{x)) = 2{n + 2£) q{x) + \x\'' Aq{x). 
Supposing further that q{x) is harmonic, we obtain that 

(3.8) A{\xWix))=2{n + 2i)q{x), 

and more generally Adxp'' q{^)) is a positive integer multiple of jxp'^"^ q{x) 
when A; is a positive integer. 

A key result now is that if p{x) is any polynomial on R", then p{x) can 
be written as a finite linear combination of polynomials of the form \x\'^'' h{x), 
where A; is a nonnegative integer and h{x) is a harmonic polynomial. To be 
more precise, if p{x) is a homogeneous polynomial of degree d, then p{x) can be 
expressed as a finite linear combination of polynomials of the form \x\'^^ h{x), 
where fc is a nonnegative integer and h{x) is a homogeneous harmonic polyno- 
mial of degree £. This can be verified using the computations above. 
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3.2 Poisson kernels, 2 

Fix a positive integer n, and let S""^ denote the unit sphere in R", so that 

(3.9) S"-i = {a; G R" : = 1}. 

Also let f{x) be a continuous real- valued function on S"~-^. We are interested 
in continuous real-valued functions h{x) on the closed unit ball -6(0, 1) which 
agree with f{x) on the boundary S"~^ and which are harmonic on the open 
unit ball B(0,1). _ 

Suppose that h(x) is a continuous real-valued function on -6(0, 1) which is 
harmonic on -8(0, 1). By standard results from advanced calculus, h{x) attains 
its maximum and minimum on -8(0, 1), and we have seen that if the maximum 
or minimum are attained in the interior -8(0,1), then h{x) is constant. As a 
consequence, the maximimi and minimum are always attained on the boundary 
S"~^, and it follows that 

(3.10) snp{\h{x)\ : x £ B{Q, 1)} = sup{|/i(a;)| : x e S"-^}, 

and in particular that h{x) = for all x e -8(0, 1) when h{x) = for all 

Thus a continuous real- valued function h{x) on -8(0, 1) which is harmonic on 
-8(0, 1) is uniquely determined by its boundary values, i.e., if hi{x), h2{x) are 
two continuous real-valued functions on -8(0, 1) which arc harmonic on -8(0, 1) 
and which arc equal on S"^^, then hi{x) — h2{x) is a continuous real- valued 
function on -8(0, 1) which is harmonic on -8(0, 1) and equal to on S"^^, and 
hence is equal to on all of -8(0,1). Let us now consider the question of 
existence of continuous functions on -8(0, 1) which arc harmonic on -8(0, 1) and 
whose values on the boundary S""'^ are prescribed in advance. When n = 1 the 
existence result is immediate, because S"~^ consists of the two points 1, — 1, 
and afRne functions are harmonic. 

In general, if f{x) is a continuous real- valued function on S"~^, then there 
is a sequence of polynomials {Pj{x)}'j^-^ on R" such that the restrictions of the 
Pj{xys to S"~^ converge uniformly to f{x). This is a well-known result from 
advanced calculus. By the results of the previous subsection, there are harmonic 
polynomials hj{x) on R" which agree with the pj(a;)'s on S""-^. 

For positive integers j. A; we have that 

(3.11) sup{|/ij(x) - /ifc(x)| : X G fi(0, 1)} 

= snp{\hj{x) - hk{x)\ : X G S""^} 
= sup{|pj(a;) -Pk{x)\ : x G S""^}. 

Using this one can check that {hj{x)}jti converges uniformly to a continuous 

function h{x) on -8(0,1) which agrees with f{x) on the boundary S""^. Also, 
h{x) is harmonic on .8(0, 1), because it follows from uniform convergence that 
h{x) is weakly harmonic. 
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One can also look for an integral formula for a harmonic extension on the 
closed unit ball of a continuous function on the unit sphere. Namely, if f{x) is 
a continuous function on S"~^, one would like to write 

(3.12) h{x)= [ P{x,y)f{y)dy 

for a harmonic function h{x) on B{0, 1) with boundary values f{x). Here dy 
denotes the usnal ckmicnt of surface integration in R". 

More precisely, one would like P{x,y) to be a continuous function defined 
for X G B{0, 1) and y G S"~^. In order for h{x) to be harmonic, P{x,y) should 
be harmonic as a fmiction of a; G B{0, 1) for each y G S"~^. In order for the 
boundary values of h to be given by / on S"~^, for each z G S"^^ the function 
P{x, y), as a function of y on S"~^, should converge to the Dirac delta function 
6z{y) at z as X G i?(0, 1) tends to z, in a suitable sense, which means in particular 
that P{x, y) should tend to for y 7^ z as a; tends to z. 

Because h{x) should be equal to 1 for all x G .8(0, 1) when f{y) = 1 for all 
y G S"~^, we should have 

(3.13) / P{x,y)dy = l 

for all X G B{0,1). If f{y) is a nonncgativc real-valued function on S""'^, 
then h{x) should be a nonnegative real- valued function on B{0,1), and this 
leads to the condition that P{x,y) should be a nonnegative real number for 
all x G B{0, 1) and y G S"-i. Also, P(0, ?y) should be equal to l/un-i for all 
y G S"~^, where denotes the area of the unit sphere S"~^. 

In fact there is a function P{x, y) with these properties, called the Poisson 
kernel for the unit ball. Specifically, P{x, y) is given by the formula 

(3.14) P{x,y)= 



The Poisson kernel is unique, and it provides another approach to the existence 
of continuous functions on the closed unit ball which are harmonic on the open 
unit ball and have prescribed boundary values. 

Similarly, on any closed ball B{w, r) a continuous function on the boundary 
sphere can be extended to a unique continuous function on the closed ball which 
is harmonic on the corresponding open ball, and the extension on the open ball 
can be expressed by an analogous Poisson integral of the boundary values. Let 
us note that the Poisson kernel is real-analytic as a function of x, which is to 
say that it admits local representations in terms of power series. Each harmonic 
function on a nonempty subset of R" is locally given as a Poisson integral, and 
as a result it follows that each harmonic function is also real analytic. 

Now let us describe an application of the existence of harmonic functions on 
balls in R" with prescribed continuous boundary values. Let U he a nonempty 
open subset of R", and let h{x) be a continuous function on U. We say that 
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h{x) satisfies the weak mean value property if for each x £ U there is a positive 
real number r{x) such that the usual mean value property holds for h for balls 
centered at x and with radius less than r{x), which is to say that 

(3.15) averageg(^_^)/i = 

for all positive real numbers r such that r < r{x). 

Suppose that C/ is a nonempty connected open subset of R", and that h{x) 
is a continuous real-valued function on U which satisfies the weak mean value 
function. If p is a point in U at which h attains its maximum, so that 

(3.16) h{x) < h{p) 

for all X £ U, then one can show that h is constant on U, with h{x) = h{p) 
for all x € U. Namely, {x E U : h{x) < h{p)} is an open subset of U simply 
because h is continuous, while {x G U : h{x) = h{p)} is an open subset of U 
because of the weak mean value property and the assumption that h attains its 
maximum at p, and hence {a; € f/ : h{x) = h{p)} = U by the connectedness of 
U. 

Another feature of the weak mean value property is that it is preserved 
by restricting a function to a smaller open set. As a result of this and the 
maximum principle in the preceding paragraph, it follows that a function which 
satisfies the weak mean value property on a nonempty open subset of R" can 
be represented locally as a harmonic function on a ball, namely a harmonic 
function whose boundary values are the given function. In short, a continuous 
function on a nonempty open subset of R" which satisfies the weak mean value 
property is actually harmonic. 

3.3 Poisson kernels, 3 

Fix a positive integer n again, and let us identify R"+^ with the Cartesian 
product R" X R, and write {x,t), x G R", i £ R, for an element of R". The 
upper half-space in R"+^ is denoted and defined by 

(3.17) Z^"+^ = {{x, i) G R" X R : i > 0}, 
and the corresponding closed half-space is given by 

(3.18) TT^^ = {(x, i) e R" X R : i > 0}. 

We shall be interested in continuous functions on the closed half-space which 
are harmonic on the open half-space. 

Let us begin with a version of the Schwarz rejection principle. Namely, if 
h{x, t) is a continuous function on the closed upper half-space which is harmonic 
on the open upper half-space, and if 

(3.19) /i(a;,0)=0 
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for all X G R", then we can extend h{x, t) to all of R" x R by setting 

(3.20) h{x,-t) = -h{x,t) 

when t > 0, and this extension is harmonic. Indeed, it is easy to see that this 

extension satisfies the weak mean value property on R" x R. 

Next, suppose that f{x) is a harmonic function on R™ which is bounded, so 
that there is a positive real number C so that 

(3.21) \.f{^)\<C 

for all X G R™. In this case one can show that f{x) is a constant function. 
Indeed, for any x,y G R™ and any positive real number r, we have that 

(3.22) f{x) - f{y) = averagc^(^_^^/ - averagc^^^. 

and one can show that the right side tends to as r ^ cx) because / is bounded. 

In fact, there is a generalization of this for harmonic functions f{x) on R"* 
which are of at most polynomial growth, in the sense that there is a positive 
real number C and a positive integer £ such that 

(3.23) \f{x)\<C{l + \x\f 

for all X G R™. In this case one can show that f{x) is a polynomial of degree 
less than or equal to For instance, one can show that the ^ih. order derivatives 
of f{x) are bounded harmonic functions on R™, and hence are constant. 

At any rate, we conclude that if hi(x,t), h2{x,t) are boimded continuous 
functions on the closed upper half-space in R™ x R which are harmonic on the 
open half-space and which are equal on the boundary, then they are equal on 
the whole closed half-space. Indeed, imder these conditions hi{x, t) — h2(x, t) is 
a bounded continuous function on the closed half-space which is harmonic on 
the open half-space and equal to when t = 0, and one can use reflection to 
get a bounded harmonic function on all of R" x R, which is therefore constant, 
and hence equal to 0. Notice that h{x, t) = t is a, harmonic function on R" x R 
which is equal to when t = and which is not equal to when t ^0. 

Now let us consider the corresponding question of existence of bounded con- 
tinuous functions on the closed half-space which are harmonic on the open 
half-space and which have prescribed boundary values. 

We can begin with the point of view of Fourier analysis. Namely, suppose 
that ^ is an element of R", and consider the bounded continuous function 

(3.24) exp(i^ • x) 
on R", as a function of x. One can check that 

(3.25) exp(- \^\t) exp(i^ ■ x) 

is a harmonic function of {x, t) on R" x R which has modulus less than or equal 
to 1 when t>0 and which agrees with the specified function of x when t = 0. 
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As in the case of the unit ball, there is a nice integral formula in this case. 
Namely, if f{x) is a bounded continuous function on R", then 

(3.26) h{x,t)= [ P{x-y,t)f{y)dy 

defines a bounded harmonic function on the upper-half space whose boundary 
values are equal to f{x), where P{x,t) is the Poisson kernel for the upper 
half-space, which is the function on the upper half-space given by 

t 

(|a;|2 + i2)(n+i)/2' 

with a„ a positive real number. 

To be more precise, if f{x) is a bounded continuous function on R", and 
we define h{x,t) on the upper half-space using the formula above, and set 
h{x,0) = f{x), then h{x,t) is a bounded continuous function on the closed 
half-space which is harmonic on the open half-space. Harmonicity of h{x, t) is 
a consequence of harmonicity of P{x, t), as a fimction on the upper half-space. 
As a function of x, P{x, t) tends to the Direc delta function at as t ^ in 
a suitable sense, and this corresponds to the fact that h{w, t) tends to f{w) as 
t 0. 

As in the case of the unit ball, P{x, t) is positive everywhere, and this corre- 
sponds to the fact that h{x, t) is positive everywhere when /(x) is nonnegative 
and not identically equal to 0. Also, 

(3.28) / P{x,t)dx = l 

JR" 

for all t>Q, which corresponds to the fact that h{x, t) = 1 for all (x, t) in the 
upper half-space when f{x) = 1 for all x G R". The Poisson kernel P{x,t) is 
the unique function on the upper half-space which represents bounded harmonic 
functions on the upper half-space that are continuous up to the boundary in 
this manner. 

3.4 Fundamental solutions 

Fix a positive integer n. Let us first consider the question of finding a Newtonian 
kernel Nn{x) in dimension n, which should be a function on R" whose Laplacian 
is equal to the Direc delta function at 0. This is to be interpreted in a weak 
sense, which is to say that 

(3.29) / Nn{x) A(l){x) dx = (piO) 

when <l){x) is a twice continuously-differentiable function on R" with restricted 
support. 

On R"\{0}, Nn{x) should be a harmonic function. This does not mean that 
Nn{x) is necessarily continuous at 0, however. It turns out that the singularity of 
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(3.27) 



Nn{x) at the origin is mild enough so that the above integral makes sense as an 
principal value integral, and behaves reasonably well, with absolute convergence 
in particular. 

The Newtonian kernel is not uniquely determined by the condition that its 
Laplacian be equal to the Dirac delta function at in the weak sense, because 
one could add to Nn{x) any harmonic function on R" and get a function with 
the same property. However, if one imposes the condition that Nn{x) be a radial 
function, then it is unique except for adding a constant. This is because radial 
harmonic functions on R"\{0} are characterized by a homogeneous second-order 
ordinary differential equation, with a two-dimensional space of solutions, and 
constant functions account for a one-dimensional space of solutions. 

When n = 1, one can take 

(3.30) N,{x) = ^\x\. 
When n = 2 one can take 

(3.31) N2{x) = b2log\x\, 
and when n > 3 one can take 

(3.32) Nn{x) = 

where 6„, n > 2, are positive real numbers. This is not too difficult to check. 

If y is a point in R", then Nn{x — y) has Laplacian, as a function of x, equal 
to the Dirac delta function at y in the weak sense, so that 

(3.33) / Nn{x- y) Acj>{x) dx = cl>{y) 

for all twice continuously-differentiable functions (j) on R" with restricted sup- 
port. Another way to look at this is that if ■i/' is a continuous function on R" 
with restricted support, then 

(3.34) *(ar)= / N^{x - y)ilj{y) dy 

is a continuous function on R" such that 

(3.35) A* = V 
in the weak sense, which is to say that 

(3.36) / ^{x)A(P{x)dx= [ il^{x) <l){x) dx 

^R" Jr" 

for all twice continuously-differentiable functions on R" with restriced sup- 
port. In general, ^I* is not quite twice continuously-differentiable when ip is 
continuous, but this is almost true, and ^' is twice continuously-differentiable 
under additional mild conditions on ij). 
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Of course the idea of a fundamental solution like this makes sense for differ- 
ential operators in general. 

For instance, when n = 1, one can consider the derivative of Ni{x), which is 
equal to —1/2 when .t < and to 1/2 when x > 0. One can make the convention 
of setting this function to be at a; = 0. At any rate, N{{x) is a fundamental 
solution for d/dx in the sense that its derivative is equal to the Dirac delta 
function at in the weak sense, so that 

(3.37) -/ N[{x)^'{x)dx = (l){0) 

Jr 

for all continuously-differentiable functions on the real line with restricted sup- 
port. 

Similarly, a nonzero multiple of 1/z defines a fundamental solution for the 
differential operator d/9z. In other words, the d/9z derivative of 1/z is equal 
to a nonzero constant times the Dirac delta function at on C, in the sense 
that 

(3.38) / -^^{z)dz 

Jc z dz 

is equal to a nonzero constant times ^(0) for all continuously-differentiable func- 
tions on C with restricted support. Note that 1/z is indeed holomorphic on 
C\{0}, as it should be to have d/dz derivative equal to there, and that I/2; is 
a constant multiple of d/dzN2{z), identifying with C in the usual manner. 
For general n, the Clifford-valued function 

En 



(3.39) En{x) = 

is a nonzero constant multiple of a fundamental solution for both and Vr. 
This means that 

(3.40) / VR(t){x)En{x)dx 
and 

(3.41) / Er,{x)VL(p{x)dx 

are equal to a constant multiple of ^{0) whenever ^{x) is a continuously- 
differentiable Clifford- valued function on R" with restricted support. This func- 
tion En{x) is a constant multiple of V^Nnix) = VnNnix), and is both left and 
right Clifford holomorphic on R"\{0}, as it should be. 

4 Harmonic and holomorphic functions, 3 
4.1 Subhcirmonic functions, 1 

In dimension 1 , the notion of subharmonicity reduces to that of convexity. Let 
(a, b) be an open interval in the real line, which is to say the set of real numbers 
X such that 

(4.1) a <x <b, 
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where of course wc assume that a < b, and also we aUow a to be — oo or a finite 
real number, and 6 to be a finite real number or +00, so that our interval may 
be unbounded, like the set of positive real numbers, or the whole real line itself. 
A real- valued function f{x) is said to be convex if 

(4.2) /(A x + {l-X)y)<X f{x) + (1 - A) f{y) 

for all X, y G (a, h) and all A e (0, 1). 

It is a well-known exercise that convex functions arc continuous. In the 
other direction, if one assumes that f{x) is a continuous real- valued function on 
(a, h), then more restrictive convexity conditions imply the general one, e.g., it 
is enough to take A = 1/2. If f{x) is a twice diffcrcntiable function on (a, 6), 
then f{x) is convex if and only if f"{x) > for all x G (a, b). 

For any continuous real- valued function f{x) on (a, b), we can say that /" > 
on (a, b) in the weak sense if 

(4.3) ( f{x)(t>"{x)dx>0 

J a 

for all twice continuously-differcntiablc functions (^{x) with restricted support 
in (a, b) such that 0(a;) > for all x S (a, b). If f{x) is also twice continuously- 
differentiable, then /" > on (0,6) in the weak sense if and only if f"{x) > 
for all X £ {a,b), basically by integrating by parts twice. For any continuous 
real-valued function f{x) on (a, 6), one can show that f{x) is convex if and only 
if /" > on (a, b) in the weak sense. 



4.2 Subhcirmonic functions, 2 

Let [/ be a nonempty open subset of R", and let f{x) be a real- valued continuous 
function on U . We say that f{x) satisifies the weak sub-mean value property 
on U if for each x gU there is a positive real number r{x) such that 

(4.4) average^^^ > f{x) 

for all r e (0,r(x)). To be more precise, this is a sub-mean value property 
for averages over balls, and we shall also consider versions with averages over 
spheres. 

By definition, if a continuous real-valued function satisfies the weak sub- 
mean value property on an open set, it also satisfies this property when restricted 
to any open subset of this set. If ?7 is a nonempty connected open subset of 
R", f{x) is a real- valued continuous function on U which satisfies the weak 
sub- mean value property, and p is an element of U at which f{x) attains its 
maximum, so that 

(4.5) fix) < f{p), 

then f{x) is constant on U. As usual, this is because {x G U : f{x) = f{p)} is 
then an open subset of U, by the weak sub-mean value property, while {x gU : 
f{x) < f{p)} is an open subset of U by continuity. 
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Let [/ be a nonempty open subset of R", let f{x) be a real- valued continuous 
function on U, let z be an element of U, and let r be a positive real number 
such that 

(4.6) r < dist(2,R"\f/). 

Let h{x) be the continuous real- valued fimction on B{z,r) which is equal to 
f{x) when jx — = r and which is harmonic on B{z,r). Thus f{x) — h{x) 
satisfies the weak sub-mean value property on B{z, r), and since f{x) — h{x) is 
equal to on the boundary, it follows that f{x) — h{x) < on B{z, r), which is 
to say that 

(4.7) f{x) < h{x) 

when |a; — 2;| < r. 

Because h{x) is harmonic on B{z,r), we have that 



(4.8) / h{x)b{x)dx = hiz) 

J B{z,r) 

whenever b{x) is a continuous real-valued function with restricted support in 
B{z,r) which is radial around z, in the sense that it depends only on |a; — 2;|, 
and which satisfies 

(4.9) / b{x) dx = 1. 

J B{z,r) 

As a result, 

(4.10) averageas(^_t)/i = h{z) 

when Q < t < r, where dB{z,t) = {x G R" : |a; — 2:| = t}. This equation also 
holds when t = r, since h{x) is continuous on B{z,r). 

It follows that 
(4-11) averageaB(^^^)/ > averageas(^_t)/ 

when <t <r, and that 

(4.12) averagers (^ ,.)/ > f{z). 

This works for all radii r < dist(z, R"\J7), and by integrating over r one gets 
that if /3(x) is a nonnegative real-valued continuous function with restricted 
support in U which is radial about z, so that (3{x) only depends on |a; — ^;|, and 
which satisfies 

(4.13) / (3{x)dx = 1, 

Ju 

then 

(4.14) / f{x)P{x)dx>f{z). 

Ju 

Also, 

(4-15) averageg^^ > f{z) 

when 2; e [/ and < r < dist(z, R"\[/). 
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4.3 Subhcirmonic functions, 3 

Let [/ be a nonempty open subset of R". A real- valued continuous function f{x) 
is said to be subharmonic if it satisfies the weak sub-mean value property on 
U, which is then equivalent to seemingly-stronger conditions, as in the previous 
subsection. If f{x) is twice continuously-differentiable on U, then Af{x) > 
for all X GU in this case. 

Conversely, suppose that f{x) is a real- valued function on a nonempty open 
subset U of R" such that / is twice continuously-differentiable and A/(x) > 
0. Actually, we can weaken this a bit by assuming that /(x) is a real-valued 
continuous function on U such that 

(4.16) ^^^0^ ^"^ (^^'^®'^^g®B(x,r)/ ~ /(^)^ - 

for all a; e ?7, or even 

(4.17) limsup r"^ ^averageg^^^^^/ - f{x)j > 

for all X gU. In this case one can check that 

(4.18) Mx)^f{x)+e\x\' 

satisfies the weak sub-mean value property on U for all e > 0. Hence 
(4-19) averageg(^_^)/e > fe{x) 

when X & U and < r < dist(a;, R"\f/), and therefore 

(4.20) average^^^^)/ > f(x) 

for the same x, r, so that f{x) is a subharmonic function on U. 

A real- valued continuous function f{x) on a nonempty open subset U of R" 
is said to be subharmonic in the weak sense if 

(4.21) / f{x) A(j){x) dx>0 

Ju 

for all nonnegative real-valued twice continuously-differentiable functions 4'{x) 
with restricted support in U. If f{x) is twice continuously-differentiable, then 
it is easy to check that this is equivalent to Af{x) > for all x & U, hy 
integration by parts. It is not too difficult to extend this argument to show that 
a subharmonic function is subharmonic in the weak sense. 

Conversely, a real-valued continuous function /(x) on U which is subhar- 
monic in the weak sense is subharmonic. This is analogous to the fact that a 
function which is harmonic in the weak sense satisfies the mean value property. 
As before, one can express differences of local radial averages of a function / 
as integrals of / times for twice continuously-differentiable functions with 
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restricted support in U, and now one should be a bit more careful to have ^ > 
too. 

Of course a real- valued continuous function f{x) on an open subset U of 

R" is harmonic if and only if f{x) and —f{x) arc subharmonic. If /i, /2 arc 
subharmonic functions on U and ai, a2 arc nonnegative real numbers, then 

(4.22) ai /i + a2 

is also a subharmonic function on U . Moreover, 

(4.23) max(/i,/2) 
is a subharmonic function on U in this case. 

4.4 Holomorphic functions and subharmonicity 

Let f7 be a nonempty open subset of R", and let f{x) be a real- valued continuous 
function on U. If a{t) is a real-valued convex function on R, then 

(4.24) a^averageg(^^)/^ < averageg^^^^a o / 

for all X & U and < r < dist(a;, R"\i7). In other words, the value of a convex 
function at an average of some numbers is less than or equal to the average of 
the values of the convex function at the same numbers. 

As a consequence, if f{x) is a real- valued harmonic function on U and a{t) 
is a convex function on the real line, then the composition a{f{x)) is a real- 
valued continuous function on U which is subharmonic. Moreover, if f{x) is 
a real- valued continuous subharmonic function on U and a(t) is a monotone 
increasing convex function on the real line, then the composition a{f{x)) is a 
real-valued continuous subharmonic function on U. For instance, in the first 
case one can take a{t) = \t\^ and in the second case one can take a{t) = when 
t < 0, a{t) = when i > 0, where p is a real number such that p>l. 

In fact, if ?7 is a nonempty open subset of C, and f{z) is a complex-valued 
holomorphic function on U, then |/(2)|p is subharmonic for all p > 0. One also 
has that log|/(2;)| is subharmonic, if one extends the notion of subharmonicity 
to allow for functions which take the value — oo, or alternatively the real- valued 
function max(log |/(2;)|, A) is subharmonic for every real number A. Indeed, 
these assertions can be verified using the fact that log |/(2;)| is harmonic on the 
set where f{z) 0. 

Now suppose that U is a nonempty open subset of R", and that f{x) is 
a C(n)-valued function on U. We can define |/(a;)| by identifying C{n) with a 
Euclidean space, using the standard basis for the Clifford algebra. It turns out 
that if f{x) is Clifford holomorphic, then |/(x)|^ is subharmonic for a range of 
p's which includes some p's strictly less than 1. There are general results of this 
type for first-order constant coefficient systems of partial differential equations 
which are "elliptic" in the sense that their homogeneous solutions are harmonic, 
and there are various situations where this applies with various ranges of p. See 
[4, 5]. 
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5 Miscellaneous 



5.1 Linear fractional transformations 

A polynomial on C is a function P{z) of the form 

(5.1) a„2" + a„_i^"-^ + --- + ao, 

where n is a nonnegative integer and Oq, . . . , a„ are complex numbers. If n > 1 
and a„ 7^ 0, or if n = 0, then we say that P{z) has degree n. If P{z) is a 
polynomial of degree n > 1, then every complex number is attained as a value 
of P{z) exactly n times if one counts multiplicities in an appropriate manner. 

Let us write C for the Riemann sphere, which consists of the complex num- 
bers C together with an extra "point at infinity" denoted 00. This is topolog- 
ically equivalent to the usual 2-dimensional sphere in the obvious manner. We 
also make the conventions that l/oo = 0, 1/0 = 00, o + 00 = 00 when a G C, 
a • 00 = 00 when a S C, a 7^ 0, etc. 

If P{z) is a polynomial on C, then P can be extended to a continuous 
mapping from the Riemann sphere in a natural way, namely, P{oo) = 00 when 
P{z) has positive degree, and P takes the same value at 00 as at other points 
when P is constant. More generally we can consider rational functions of the 
form 

(5.2) "^"^'wr 

where Q{z) is a polynomial which is not identically equal to 0. To be more 
precise, R{z) is then defined as a complex number when z is a complex number 
and Q{z) ^ 0, while if Q{z) = or 2; = 00 one can follow the usual conventions, 
so that the rational function remains continuous as a mapping from the Riemann 
sphere to itself. 

Assuming that the polynomials P{z), Q{z) have no common factors, then the 
degree of the corresponding rational function R{z) is defined to be the maximum 
of the degrees of P{z) and Q{z). If R{z) is a rational function of degree n > 1, 
then each element of the Riemann sphere occurs as a value of R{z) n times, 
with appropriate multiplicities. This is not too difficult to see. 

A linear fractional transformation on the Riemann sphere is a rational func- 
tion of degree 1. Explicitly, this is a mapping which can be expressed as 

(5.3) 

cz + d 

where a, b, c, d are complex numbers such that 
(5.4) ad-bcj^O. 

This condition ensures that the denominator in the linear fractional transfor- 
mation is not identically equal to 0, and that the numerator and denominator 
are not simply constant multiples of each other. Notice that if we multiply a, b, 
c, and d by the same nonzero complex number, then this admissibility condition 
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is still satisfied, and that the corresponding linear fractional transformation on 
the Riemann sphere is in fact the same as the original one. 

In fact, we can think of the complex numbers a, b, c, d as entries of a 2 x 2 
matrix 

'a b' 
c d. 



(5.5) 



whose determinant is then equal to a c — 6 d, so that the admissibility condition 
is equivalent to the requirement that this matrix be invertible. One can check 
that the composition of two linear fractional transformations is again a linear 
fractional transformations, and that this corresponds to multiplying the associ- 
ated matrices in the usual way. In particular, the inverse of a linear fractional 
transformation can be given by the inverse of the associated matrix. 

It is helpful to consider some special cases of linear fractional transforma- 
tions, of the form 

(5.6) az 
for a nonzero complex number a, 

(5.7) z + b 

for any complex number b, and 
(5.8) 
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One can check that every linear fractional transformation can be expressed as 
a combination of these basic types. Of course this can also be viewed in terms 
of 2 X 2 matrices. 

Let us consider this in another way, which extends to higher dimensions. 
Namely, if n is a positive integer, let CP" denote the space of complex lines 
through the origin in C""'""'^, the space of (n + l)-tuplcs of complex numbers. 
In other words, we can start with C"+^\{0}, and identify two elements z, w of 
this space when there is a nonzero complex number A such that w = Xz. 

If A is an invertible complex linear mapping on C"+-'^, then A takes the 
origin in C""*"^ to itself, and it takes complex lines in C""*"^ to themselves, 
and thus induces a transformation on CP". Multiplication of ^4 by a nonzero 
complex number does not change the induced transformation on CP". Compo- 
sitions of invertible linear mappings on C"^^ correspond to compositions of the 
induced transformations on CP", and in particular the inverse of an induced 
transformation on CP" is induced by the inverse of the original linear mapping 
on C"+i. 

There is a natural way to embedd C" into CP", which is to take an element 

z of C" and extend it to a nonzero clement of C"+^ by adding an {n + l)th 
coordinate which is set equal to 1. In other words, we identify C" with the 
affine hyperspace in C"+^ of points whose last coordinate is equal to 1, and 
then each point in this affine hyperspace determines a complex line in C"^"'^ 
through the origin, namely the complex line that passes through that point. 
This accounts for all the complex lines in C""*"^ except for those which lie in the 
linear subspace of points whose last coordinate is equal to 0. 
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Roughly speaking then, C" looks like C" together with a copy of CP""'^, at 
least when n > 2. When n = 1, the embedding of C into CP^ accounts for all 
the lines in through the origin except for one line, which is the set of points 
with second coordinate equal to 0. Thus CP^ looks like C with one additional 
element, which gives another way to think about the Riemann sphere. 

Linear fractional transformations on the Riemann sphere are then the same 
as transformations on CP"'^ induced by invertible complex linear transformations 
on C'^. To be more explicit, an invertible linear transformation on C'^ can be 
given as 

(5.9) {z,w) {a z + bw,cz + dw), 

where again a, 6, c, d are complex numbers such that ac — bd 0. This leads 
to the earlier formula through the identifications described above. 

There is another and rather different way to extend the idea of linear frac- 
tional transformations to higher dimensions. For this we think of as being 
identified with C in the usual way. Now we work on R", and on R" ~ R"U{(X)}, 
which can be identified topologically with an n-dimensional sphere. 

To describe this generalization, one can go back to the idea of building 
blocks. Namely, one can use translations, dilations, orthogonal transformations, 
and inversion about the unit sphere in R" as the basic building blocks. One 
then gets a nice family of transformations on R" generated by these. 

Actually, this would correspond to also allowing complex conjugation as a 
transformation on the Riemann sphere. Alternatively, in n dimensions, one can 
add the requirement that the transformations preserve orientations. At any rate, 
this generalization is closely connected to conformal geometry in n dimensions. 



5.2 Power series and the exponential function 

Consider a power series 

oo 

(5.10) ^anz", 

n=0 

where the a„'s are complex numbers and z is a complex variable. We shall be 
interested in the set of 2: G C for which this series converges. Of course this 
series converges trivially when z = 0. 

If the series converges for some Zq & C, then we have that the sequence 
{a„ Zq}^^q tends to 0, and is bounded in particular, so that there is a nonneg- 
ative real number C such that 

(5.11) \an\\zor<C 

for all n. Assume that 2^0 7^ and that z is a complex number such that 

(5.12) \z\ < \zo\. 
Under these conditions, we have that 

(5.13) \an\\zr<C 




46 



and hence Y^^^q an -z" converges absolutely by comparison with a convergent 
geometric series. 

Using this simple remark, one can check that either the power series con- 
verges only for z = 0, or there is a positive real number R such that the scries 
converges absolutely when \z\ < R and does not converge when \z\ > R, or the 
series converges absolutely for all ^ e C. In the first case we can set ii = 0, and 
in the third case we can set = oo. We call R the radius of convergence of the 
power series. 

If the power series converges absolutely for some complex number zq, then 
it also converges absolutely for all complex numbers z such that \z\ < \zo\. In 
general, if a power series has a radius of convergence R which is positive and 
finite, then the series might converge for no complex numbers z with \z\ = R, 
or for some of them, or for all of them. For instance, the series 



(5.14) 

n=0 

has radius of convergence equal to 1 and converges for no ^ € C with 1^1 = 1, 



n=l 



has radius of convergence 1 and converges for ^; e C with |z| = 1 and z / 1, 
and does not converge when z = 1, and 



oo ^ 
Z^ 



(5-16) 



n=l 

has radius of convergence 1 and converges absolutely for all ^; e C with \z\ < 1. 
Let 

OO OO 

(5.17) 

n=0 n=0 

be two series of complex numbers, and define an associated series 



(5.18) ^C7„ 



n=0 

by 

n 

(5.19) Cn = J2^j^n-j. 

This is called the Cauchy product of the original two series. Formally, the sum 
of the C„'s is equal to the product of the sums of the A„'s and the i?„'s, and 
indeed this is the case if the A„'s and Bn's are equal to for all but finitely 
many n. If the original series are power series, with A„ = a„ Zn and -B„ = 6„ z"', 
then Cn = Cn -z", with c„ = J^j^o "-J ^n-j- 



47 



It is not too difficult to show that if ^ An , converge absolutely, then 

the Cauchy product series J2 also converges absolutely, and the sum of the 
Cauchy product is equal to the products of the sums of the initial series. One 
can also show that if both scries converge and at least one converges absolutely, 
then the Cauchy product converges, and the sum of the Cauchy product series 
is equal to the product of the sums of the original series. It is not true in general 
that the Cauchy product scries converges when the two initial scries converge. 

There is a nice result though which says that if the two initial series converge 
and the Cauchy product series converges, then the sum of the Cauchy product 
series is equal to the product of the sums of the two initial scries. To show this, 
one can use the notion of Abel summability. Basically, it is convenient to view 
the series as power series. 

Specifically, for a given series J2'^=o associated Abel sums are defined 

by 

oo 

(5.20) 5]r"A„, 

where r is a positive real number such that r < 1. We assume that these 
series converge for all r E (0, 1), which is the same as saying that they converge 
absolutely for all r e (0, 1), as before. This holds under modest assumptions on 
the size of the An's. 

The series X^^q said to be Abel summable if the Abel sums exist for 
all r S (0, 1), and if the limit of the Abel sums exist as r — > 1—. This limit is 
then called the Abel sum of the series Y^'^^q A^. For instance, if the A„'s are 
nonnegative real numbers, then the Abel sums are monotone increasing as a 
function of r, and it is easy to see that the series is Abel summable if and only if 
it converges in the usual sense, and the Abel sum is equal to the ordinary sum. 

If J2 converges absolutely, then it is again easy to see that J2 is Abel 
summable, and that the Abel sum is equal to the ordinary sum. It turns out 
that this works in general when Y An converges. One can approach this by 
expressing the Abel sums in terms of the partial sums of the original series, 
noting also that the Abel sums converge for all r G (0, 1) because the ^„'s are 
bounded. 

There are also series which are Abel summable and which do not converge 
themselves. A basic family of examples is given by the series 



(5.21) 



n=0 



where ^; is a complex number such that = 1 and z ^ 1. When r G (0, 1), we 
have that 



1 

(5.22) Z1^"'^" = Y3" 



rz 

n=0 



and this converges to (1 — z)^^ as r ^ 1 — . 

Now suppose that A„, '^Bn are two convergent series of complex num- 
bers. In particular their terms are bounded, and the Cauchy product series 
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has the property that C„ = 0(n + 1), say, which is sufficient to ensure 
that X^r" C„ converges for ah r G (0, 1). If we assume that X]C„ converges, 
then the sum of this series is the same as the limit of the Abel sums r" C„, and 
these Abel sums arc equal to the product of the Abel sums ^ r" A„, ^ r" _B„, 
since we have absolute convergence when r € (0, 1), and it follows that J2 is 
equal to ^ A„ times ^ -B„. 

Suppose that we have a power series 

oo 

(5.23) ^a„^" 

n=0 

which converges for z G C with \z\ < i? for some i? > 0. A basic result states 
that if < r < -R, then the partial sums of this series converge uniformly to the 
whole sum for a; e C such that < r. In particular, it follows that the power 
series defines a continuous function on the disk 

(5.24) {zeC:\z\< R}. 
Under these conditions we also have that the series 



(5.25) X^na„ 



n=0 



converges for \z\ < R. The original series defines a holomorphic function f{z) 
on the open disk in the complex plane with center and radius R, and this 
series represents the complex derivative f'{z) of f{z). 

Now let us look at a specific power series. Namely, consider the series 



n=0 

where as usual n! denotes "n factorial" , the product of the integer from 1 to n, 
which is interpreted as being equal to 1 when n = 0. For each positive integer 
k, we have that 

(5 27) 1 <^^ ^ 

^ ' ' n! - (fc-l)!A;"-fe 

when n > k, and one can use this to show that our power series converges for 
all z in C. 

Of course this is the series for the exponential function exp(2;). With this 
definition one can see directly that 

(5.28) exp{z + w) = exp(z) exp{w) 

for all z,w S C. This is an example of the general result for Cauchy products 
of absolutely convergent series, as one can check using the binomial theorem. 

Clearly exp(O) = 1. When z is a positive real number, exp(2:) is a positive 
real number greater than 1. If 2; is a negative real number, exp(z) = 1/ exp(— 2;) 
is a positive real number less than 1. 
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For any complex number z, 



(5.29) exp(2;) = exp(^;), 

where a denotes the complex conjugate of a complex number a. In particular, 
\i z = X + iy, where x, y are real numbers, then 

(5.30) |exp(2;)| = exp(a;), 
which is to say that 

(5.31) |exp(iy)| = l. 
In fact, 

(5.32) exp{iy) = cosy + i siny, 

as one can see from the power series expansions. 

The exponential function is a holomorphic func;tion of z whose complex 
derivative is equal to itself, just as the ordinary derivative of exp(2;) for x G R 
is equal to itself, as one can see from the power series. Note that one can also 
derive relations for the derivatives of the real and imaginary parts of exp(j y) as 
a function on the real line which are compatible with the formula in terms of 
sines and cosines above. The fact that the exponential of a sum is equal to the 
product of the individual exponentials also contains standard addition formulas 
for sines and cosines. 

As a consequence of the trigonometric relations, we get that 

(5.33) exp(27rin) = l 

when n is an integer, and that these are the only times when exp(2:) = 1. Also, 
exp{z) 7^ for all z G C, and if C is a nonzero complex number, then there are 
z's in C such that 

(5.34) cxp(z) = C- 

Specifically, the real part x of z is uniquely determined by the condition that 
exp(a;) = |C|, while the imaginary part y of z should be taken to be an "angle" 
for C, as in polar coordinates for ~ C, which is determined up to adding 
integer multiples of 27r. 

If J7 is a nonempty open subset of C\{0} and A(C) is a continuous complex- 
valued function on U, then we say that A(C) is a branch of the logarithm on U 
if 

(5.35) exp(A(C)) = C 

for all ( &U. Assuming that Ai(C), A2(C) are branches of the logarithm on the 
same nonempty open subset U of C\{0}, we get that 

(5.36) Ai(C)-A2(C) e27riZ 

for all ( G U. Under the additional hypothesis that U is connected, it follows 
that there is a single integer n such that 

(5.37) Ai(C)-A2(C) = 27rin 
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As above, if X{Q is a branch of the logarithm on a nonempty open subset U 
of C\{0}, then 

(5.38) ReA(C)=log|C|, 

where the logarithm logp of a positive real number p is the real number r 
characterized by exp(r) = p. Also, A(C) is a holomorphic function on U. The 
complex derivative of X{Q is given by 

(5.39) A'(C) = J. 

The function 1/^, up to a nonzero constant multiple, is a fundamental solu- 
tion foTc d/dz on C. On R" for general n one has the Clifford- valued function 

(5.40) 1^2=1^ 

a nonzero constant multiple of which is a fundamental solution for D^, Vji. 
Integrals of this function are also related to measurements of angular increments, 
as in the complex plane. 

5.3 Linear algebra in R", continued 

Let Z7 be a nonempty open subset of R", and let f{x) be a real-valued twice 
contimiously-differentiable function on U . For each a; G ?7, we get the matrix of 
second partial derivatives of /, 

and a theorem of advanced calculus states that this matrix is symmetric. The 
Laplacian Af{x) of / at a; is the same as the trace of this matrix. 

If T is a linear transformation on R" and t; is a nonzero vector in R", then 
V is said to be an eigenvector for T with eigenvalue A G R if 

(5.42) T{v) = Xv. 

If T is a symmetric linear transformation on R", and if v is an eigenvector of 
T with eigenvalue A, then T maps every vector w in R" orthogonal to t; to a 
vector orthogonal to v. Indeed. 

(5.43) {T{w),v) = {w,T{v)) 

= X{w,v)=0. 

In general, a linear transformation T on R" is said to be diagonalizable if 
there is a basis of R" consisting of eigenvectors of T, which is equivalent to saying 
that there is an invertible linear transformation A on R" so that the matrix 
associated to AT is diagonal. If T is a symmetric linear transformation on 
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R", then there is an orthonormal basis of R" consisting of eigenvectors of T, 
which is equivalent to saying that there is an orthogonal linear transformation 
A on R" such that the matrix associated to AT is diagonal. To find an 
eigenvector v for T, one can maximize 

(5.44) {T{v),v) 

among vectors v such that \v\ = 1, and then one can repeat the process to get 
an orthonormal basis of eigenvectors. 

Suppose that T is any linear transformation on R". Then T* T is a nonneg- 
ative symmetric linear operator on R", and there is a nonnegative symmetric 
linear operator T on R" such that 

(5.45) f 2 = T. 

Namely, there is an orthonormal basis of R" consisting of eigenvectors of T* T, 

with nonnegative eigenvalues, and one can define T so that the same basis is a 
basis of eigenvectors and the eigenvalues are nonnegative square roots of those 
of T. 

Observe that 

(5.46) \f{v)\ = \T{v)\ 

for all V G R", since 

(5.47) \f{v)\' = {nv),fiv)} 

= {{T*T){v),v) = {Tiv),Tiv)) = \Tiv)\^ 

Using this one can show that there is an orthogonal linear transformation A on 
R" such that 

(5.48) T = Aof. 

This is called a polar decomposition of T, although one might prefer a decom- 
position of T in the other order, as the composition of a nonnegative symmet- 
ric operator with an orthogonal transformation, and this can be obtained by 
expressing T* as the composition of an orthogonal transformation and a non- 
negative symmetric operator, and taking the adjoint afterwards. For this polar 
decomposition, one can simply take A = T when T is invertible, so that T is 
invertible, and otherwise one has some freedom in choosing A on the orthogonal 
complement of the image of T. 

Now let us look some more at orthogonal transformations in general. If 
Vi, . . . ,Vn and wi, . . . , are orthogonal bases on R", then there is a imique 
linear transformation A on R" such that A{vj) = Wj for 1 < j < n, and A is an 
orthogonal transformation. Conversely, if A is an orthogonal linear transforma- 
tion on R" and vi, . . . , w„ is an orthonormal basis of R", then A{vi), . . . , A{vn) 
is also an orthonormal basis of R". 

If A is an orthogonal linear transformation on R" and Ui , . . . , u„ is an or- 
thonormal basis of R", then A is almost determined by its values onvi,. . . , Vn-i- 
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Namely, A{vn) should be a unit vector orthogonal to A{vi), . . . , and 
there are only two choices for this vector, which are negatives of each other. If 
^ is a rotation, which is to say an orthogonal linear transformation with deter- 
minant 1, so that A preserves orientations on R", then A is determined by its 
values on vi,. . . , Vn-i- 

Suppose that u\,. . . ,Uk are orthonormal vectors in R" with fc < n — 2, or 
simply that n > 2 if no such vectors are specified. If v, w arc unit vectors 
in R" which are orthogonal to ui,...,Uk, then one can make a continuous 
deformation of t; to «; through unit vectors in R" which are orthogonal to 
Ml, . . . Indeed, there is a continuous deformation of the identity operator 
on R" through orthogonal linear transformations on R" which fixes Uj for 
j = 1, . . . ,k and which deforms v to w. 

Using this, one can show that any two rotations on R" can be continuously 
deformed to each other through orthogonal linear transformations. For that 
matter, any two invertible linear transformations on R" whose determinants 
have the same sign can be continuously deformed to each other. If Ai , A2 are two 
invertible linear transformations on R" whose determinants have opposite signs, 
then any continuous deformation from Ai to A2 through linear transformations 
on R" must pass through at least one linear transformation with determinant. 

Recall that a function f{x) on R" is said to be radial if it depends only on 
This is equivalent to saying that 



for all orthogonal linear transformations A on R". For ii x, y are elements of 
R" such that |x| = then there is an orthogonal linear transformation A on 
R" such that A{x) = y. 

Suppose that C/ is a nonempty open subset of R" and that / is a twice 
continuously-differentiable function on U. If A is an orthogonal linear trans- 
formation on R", then f o A is a twice continuously-differentiable function on 
A~^{U). One can check that the Laplacian of / o A at a point x G A~^{U) is 
equal to the Laplacian of / at A{x). 

In particular, / o ^ is harmonic on A^^{U) if / is harmonic on U. On the 
complex plane, if </> is a holomorphic function on an open set V which takes 
values in another open subset U of C, and if / is a harmonic function on U, 
then / o (/) is a harmonic function on V. If / is also holomorphic, then f o cj) is 
holomorphic on V. If /(z) is a holomorphic function on an open subset U of C, 
then 



is a holomorphic function on ?7 = {z : z G U}. 

Let A be an orthogonal linear transformation on R" with associated matrix 
(aj^k), so that the jth coordinate of A{x) can be written as 



(5.49) 



f{A{x)) = fix) 



(5.50) 



n 



(5.51) 
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Wc can associate to A a linear isomorphism A on the Clifford algebra C{n), 
defined by the conditions that A(r) = r when r is a real number, 

n 

(5.52) A{ek) = y^,aj,kej, 

and of course A acting on a product of e;'s is equal to the product of the 

corresponding A(e;)'s. One can check that the assumption that A be an orthog- 
onal transformation provides the appropriate compatibility conditions, namely, 
A{ei)'^ = — 1 for / = 1, . . . , n and A{ei) A{em) = —A{em) (e;) when I ^ m, which 
basically means that these products have real part equal to 0. 

Suppose that U is a. nonempty open subset of R", and that f{x) is a C{n)- 
valued continuously differentiable function on U, and let /a denote the function 
on A{U) defined by 

(5.53) fAix) ^ A{fiA-\x))). 

One can check that Vl or applied to Ja at a point x G A{U) is the same as 
A applied to or Vn of /, respectively, evaluated at A~^{x). In particular, 
Ja is left or right Clifford-holomorphic on A{U) if / is on U. 

It is easy to see from the definition of Clifford-holomorphicity that at each 
point in the domain of a Clifford-holomorphic function, the derivative of the 
function in the direction of one of the coordinates is determined by the deriva- 
tives of the function in the directions of the other coordinates. Using the invari- 
ance under orthogonal transformations described in the previous paragraph, one 
can check that the derivative of a Clifford-holomorphic function in any direction 
is determined by the derivatives in the orthogonal directions. Of course this is 
very clear in the classical case of holomorphic functions of a single complex 
variable. 

We can also derive as follows. Let u be a vector in R", v = X^jLi Vj ej the 
corresponding element of C{n), and f{x) a C(n)-valued continuously differen- 
tiable function on an open subset U of R". If VLfix) = for some x gU, then 
of course 

n p. 

(5.54) y^ei—fix)=0, 

which can be rearranged as an identity between the directional derivative of / 

at X in the direction of v and a combination of terms of the form Vj d/dxif{x) — 
vi d/dxjf{x) times ej e;. These terms involve directional derivatives of / at a; in 
directions orthogonal to v, and an analogous computation works for 'Dfif{x)v 
when Diifix) = 0. 

5.4 Some loose ends 

Let U he a nonempty open subset of R". If /i, /2 are two real- valued continu- 
ously differentiable functions with restricted support in U, consider the bilinear 
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form S(/i,/2) acting on /i, /2 defined by 

(5.55) B{h,h) = [ {yh{x),Vh{x))dx, 

where V f{x) denotes the gradient of / at x, the vector whose components are 
the partial derivatives 

(5.56) 

Of course B{fi,f2) is symmetric in /i, /2, and if /i is also twice continuously- 
differentiable, we have that 

(5.57) B{h,h) = - I Ah{x)f2{x)dx. 

Ju 

Also notice that -B(/, /) > 0, and indeed B{f, /) = if and only if / = 0, since 
we are considering functions with restricted support. 

One can extend this to consider functions which tend to at the boundary 
of U and which have enough regularity, and it can be of interest to make some 
assumptions about U as well. In particular, for a nice domain U, one can try 
to diagonalize this bilinear form, and thus the Laplacian, with respect to the 
standard integral innner product 

(5.58) / h{x)h{x)dx, 

Ju 

and for functions which tend to at the boundary of U . This means looking 
for smooth functions f on U which are eigenfunctions of the Laplacian, so that 
there is a real number A such that 

(5.59) Af{x) = Xf{x) 

for all x gU, and which tend to at the boundary of U. 

Just as for harmonic functions, one can define what it means for a continuous 
function on U to be an eigenfunction for the Laplacian in the weak sense, and 
show that this implies that the function is smooth. There are also basic theories 
in analysis for showing that there are eigenfunctions of the Laplacian which tend 
to at the boundary, and enough of them to diagonalize the Laplacian among 
functions which tend to at the boundary. Of course there are some details in 
formulating this, and indeed some variations along these themes. 

If h{x) is a real- valued harmonic function on a nonempty open subset of R", 
and if (j){x) is a real-valued continuously difi^erentiable function with restricted 
support in U, then 

(5.60) / \V{h + (p){x)fdx= I \Vh{x)fdx+ [ \V(l){x)f dx, 
Ju Ju Ju 

since 

(5.61) / {Vh{x),V(j){x))dx = - I Ah{x) (pix) dx = 0. 
Ju Ju 



55 



In particular, notice that 



(5.62) 



Ju Ju 



[ \V{h + <l)){x)\^ dx> [ \\/h{x)\^dx. 



More precisely, this works under suitable conditions for any function (j) which 
vanishes at the boundary of U. Conversely, one can use this as a basis for 
looking for harmonic functions on U with given boundary values, by trying to 
minimize 



among functions with those boundary values. 

Another approach to finding harmonic functions on a nonempty open subset 
U of R" with specified boundary values is to consider the supremum of sub- 
harmonic functions on U which are less than or equal to the prescribed values 
at the boundary. A third method is to generate a lot of harmonic functions on 
U in general, and then choose ones with desired properties. We have already 
seen "fundamental solutions" on R"\{0} which are harmonic or holomorphic 
and hence harmonic, and one can use translates of these to get harmonic func- 
tions on R"\{j/} for each y in R", and then linear combinations of these with 
y e R"\{7 to get harmonic functions on U. 
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